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Abstract 

We construct composite CFT operators from a large number of fermionic primary 
fields corresponding to states that are holographically dual to a zero temperature Fermi 
gas in AdS space. We identify a large N regime in which the fermions behave as free 
particles. In the hydrodynamic limit the Fermi gas forms a degenerate star with a radius 
determined by the Fermi level, and a mass and angular momentum that exactly matches 
the boundary calculations. Next we consider an interacting regime, and calculate the 
effect of the gravitational back-reaction on the radius and the mass of the star using 
the Tolman-Oppenheimer-Volkoff equations. Ignoring other interactions, we determine 
the " Chandrasekhar limit" beyond which the degenerate star (presumably) undergoes 
gravitational collapse towards a black hole. This is interpreted on the boundary as a 
high density phase transition from a cold baryonic phase to a hot deconfined phase. 
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1 INTRODUCTION 

String theory contains all kinds of elementary particles and forces, and thus in principle can 
be used to describe all forms of matter. In the past two decades much of the progress in string 
theory has centered around black holes, which may be viewed as the most extreme form of 
matter. In fact, the matter that makes up a black hole is in such an extreme form that it has 
entirely disappeared from the picture and is replaced by pure geometry. Holography may be 
seen as an attempt to rescue some of this matter by assigning material properties to the black 
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hole horizon. To clarify the relationship between the holographic matter and the matter that 
made the black hole one needs to investigate the process of black hole formation, in other 
words, gravitational collapse. 

A related question is whether string theory, in the context of the AdS/CFT correspondence 
PP, can be used to describe other extreme forms of matter, such as the degenerate matter 
that makes up a neutron star. Indeed, neutron stars, or more generally degenerate stars like 
"quark stars" or "strange stars", can be seen as black hole precursors in the sense that they 
may undergo gravitational collapse to form a black hole. An essential feature of neutron stars 
is that they are made out of fermions: it is the exclusion principle that keeps the star from 
undergoing gravitational collapse. Our aim is to obtain a description of a "neutron star" (or 
more precisely a "degenerate star") in anti de Sitter space in terms of the dual conformal 
field theory on the boundary. In particular, we are interested to investigate and interpret the 
criterion that determines whether the star undergoes gravitational collapse. 

One of the first and most famous equations that deals with the physics of a neutron star 
is the Tolman-Oppenheimer-Volkoff (TOV) equation. In four dimensional flat space time it 
takes the form 

dp{r) 2 £, £>(r) + p(r) M(r) + Anr 2 p(r) 
dr ~ ~ 7 2 1 - 2GM{r) Jr ^ ' ' 

where the radial dependent mass M(r) is obtained from the energy density p(r) through 

— -±-L = Airr z p(r). (1.2) 
dr 

The TOV equation follows from combining the energy-momentum conservation with the Ein- 
stein equations and holds for spherically symmetric neutron stars. It describes the radial 
dependence of the pressure p(r) and energy density p(r) for a given equation of state. Solu- 
tions only exist when the pressure is sufficient to sustain the gravitational force to prevent 
the neutron star from collapsing. Typically one finds that there is an upper limit for the mass 
beyond which the neutron star will be too heavy and starts to collapse to form a black hole 
(or in some cases another type of degenerate star). In their original paper [2] Oppenheimer 
and Volkoff used the equation of state of a free relativistic Fermi gas to find a limiting mass 
of 0.7 times the solar mass. However, more realistic equations of state lead to a higher value 
of the limiting mass around 2 or 3 solar masses. 

In this paper we construct states in the boundary CFT that are the holographic dual of a 
degenerate star in anti de Sitter space. Our discussion will be quite general without reference 
to a specific AdS/CFT setup and the qualitative aspects of our results are independent of 
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the number of space time dimensions. Using the state-operator correspondence we consider 
degenerate states corresponding to composite multi-trace operators constructed from a large 
number of copies of a given fermionic single trace operator. Using 't Hooft factorization it is 
argued that in a strict large N limit these operators behave as (generalized) free fields. The 
bulk description of these states is given by a zero temperature degenerate gas of relativistic 
fermions. In the limit of a large number of particles one can give a hydrodynamic description 
purely in terms of their energy density and pressure. We first consider the case without bulk 
interactions so that the equation of state is simply that of a free degenerate Fermi gas. We 
show that in the hydrodynamic description the degenerate star is confined to a spherical 
region in AdS space, due to the gravitational potential well. The radius in AdS units turns 
out to be given by the ratio of the Fermi momentum in the center of the star and the mass of 
the constituent fermions. As a verification of our hydrodynamic description we show that it 
precisely reproduces the mass predicted from a boundary analysis for static degenerate stars 
as well as for the case with rotation. 

Next we consider a double scaling regime chosen so that the gravitational self interactions 
of the fermions will become important, while we can still ignore Planck scale physics. Since 
our primary interest is to understand first qualitatively the criterion for gravitational collapse 
we focus our attention to including the self gravity of the degenerate star, while ignoring the 
role of all other interactions. Thus, following Oppenheimer and Volkoff, we use the equation 
of state of a free degenerate Fermi gas in the hydrodynamic description of the star, and take 
the gravitational self interaction into account through a generalization of the TOV equation 
to five dimensional anti de Sitter space. The resulting equations are studied numerically. 
We find that as a function of the central density, the mass of the star at first increases, but 
reaches a limiting mass at some critical value after which the mass slightly reduces again. 
We interpret the limiting mass as the point after which the star will become unstable and 
will start to undergo gravitational collapse. The value of the limiting mass translates in the 
boundary theory to a limiting conformal dimension of the composite operator made from the 
fermionic primary fields. 

A basic motivation for our work was towards analyzing the process of gravitational col- 
lapse and black hole formation in AdS/CFT. This is an interesting problem, which has been 
addressed in several works in the past [3HT6]. The basic intuition is that radial collapse 
corresponds to motion in the "scale direction" of the CFT, while the subsequent black hole 
formation corresponds to deconfinement and thermalization on the boundary theory. A more 
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quantitative description of this picture would clearly be desirable, but is difficult since the 
boundary theory is strongly interacting, the process time dependent and the radial holo- 
graphic dimension is not easily represented in the CFT. In this direction, it may be helpful 
to identify simple initial conditions to study (at least the onset of) the dynamical process of 
collapse. A star, on the verge of gravitational collapse, is a good candidate since it is a static 
configuration and presumably simpler to understand from the boundary point of view. 

Another motivation was to make contact with the recent developments on the holographic 
analysis of strongly interacting fermions. In several works [TTH20] correlators of fermionic fields 
have been evaluated in the presence of a black hole in AdS. Such computations reveal the 
presence of a Fermi surface in the bulk, of non-Fermi liquid type. For a better understanding 
of these results it might be worthwhile to explore the boundary meaning of a simpler system, 
that of a Fermi-gas in AdS without the presence of a black holqj- One should keep in mind that 
while the fermions are free in the bulk, they are not ordinary free fermions on the boundary 
(i.e. they do not obey a Dirac equation on the boundary) as the boundary theory is strongly 
coupled. Instead they are "generalized free fields" i.e. their correlators factorize even though 
they do not obey linear wave equations. 

Finally let us mention another motivation for constructing a holographic degenerate star, 
one of a more conceptual nature. Our degenerate star in AdS resembles, for many purposes, 
a conventional macroscopic object. Since it is embedded in AdS it has a holographic descrip- 
tion on the boundary theory. It would be interesting to understand in more detail how the 
hologram works in this case. In particular it might be easier to "decode" the hologram in 
the case of a star than for a black hole, due to the absence of a horizon. Our star is a static, 
spherically symmetric configuration in the bulk. The only nontrivial information is encoded 
in the radial profile of the density. It is not obvious how this radial profile is encoded in 
the boundary theory, since the radial dimension is not manifest on the boundary but rather 
entangled with the "scale" of excitations in the CFT. We believe that a star in AdS is a good 
toy model to further explore the holographic mapping. 

We would like to emphasize the following point about the aforementioned goal. One 
might think that decoding the hologram on the boundary theory is difficult because of strong 
coupling effects in the gauge theory and thus hopeless unless one manages to compute at 
strong coupling. While this is true to some extent, we believe that there are interesting 
questions to be understood, which are not dependent on strong coupling dynamics. In a large 



5 The interesting paper |21j appeared while this draft was being prepared, in which fermion correlators 
were evaluated on a background with fermions, without a black hole horizon. 
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N gauge theory the strong coupling effects can be "hidden" in the conformal dimensions and 
3-point functions of single-trace gauge invariant operators. These quantities are difficult to 
compute. Nevertheless, let us assume that in some magical way we knew their exact values 
at strong coupling. Would that be sufficient to understand how the hologram works? In 
principle it shoulqj, because from this information we can reproduce all correlators of the 
CFT. However in practice, even armed with this information, we do not know yet how to 
reconstruct the bulk. Understanding how the radial profile of a macroscopic object, such as 
our degenerate star, is encoded on the boundary theory may be a good starting point for 
addressing this question. Ultimately it would be fascinating to understand the CFT meaning 
of the Einstein/TOV equations^. Unfortunately we do not have much to report on this yet 
but we hope to revisit this and related questions in future work. 

The main idea presented in this work was briefly described in [22]. In this paper we 
expand on various aspects that were not covered in detail in [22], we consider the case of 
charged fermions and that of a thermal star, and finally we discuss some issues about the 
validity of our approximations and the embedding in specific AdS/CFT dualities. The plan 
of the paper is as follows: in section |2] we consider a ball of Fermi gas in AdS in the limit 
of infinite N and construct its dual holographic representation. In section [3] we discuss our 
double scaling limit, estimate the importance of various interactions and present the TOV 
equations for the fermionic fluid. In section H] we present our numerical results and discuss the 
Chandrasekhar limit beyond which the star undergoes gravitational collapse. In section [5] we 
try to understand the boundary interpretation of the gravitational interactions between the 
fermions. In section |6] we discuss the extent to which our various approximations are justified 
and the possibility of embedding our construction in a specific AdS/CFT correspondence. In 
section [7] we analyze the case where the fermions are charged under a U(l) gauge field in the 
bulk. In section [8] we present a similar system corresponding to a star made out of thermal 
gas and discuss its relevance as the dual of a "superheated phase" in certain gauge theories. 
Finally in section [9] we close with some discussions. 



6 At least for length scales which are not too small compared to the size of AdS and in the absence of black 
holes. In other words when the 0(N 2 ) degrees of freedom are not excited nontrivially above their ground 
state at strong coupling. 

7 In CFTs with classical gravity dual. 
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2 FREE FERMI GAS IN ADS 



In this section we will study a degenerate gas of free fermions in AdS^+i- We will assume that 
the bulk theory is holographically dual to a conformal field theory which has the analogue of 
a large N expansion. To be more general we will express the expansion parameter in terms 
of the central charge c of the CFT. In this section we work in the limit c — > oo, keeping other 
quantities fixed. In this limit all interactions between single-trace operators are suppressed by 
powers of 1/a/c and the fermions become essentially free. The only interaction which remains 
is the gravitational attraction of the AdS background which acts as a confining box. We work 
in a regime of the gauge theory where there is a classical gravity dual, such as the A 3> 1 
limit in the J\f = 4 SYI#. 

2.1 Boundary description 

We will construct composite multi-trace operators that are made out of a large number of 
single trace fermionic operators. The composite operators will be "degenerate" in the sense 
that they represent the operators with the lowest possible conformal dimension that can be 
made out of a given number of single trace fermionic operators of a certain kind. Through 
the state-operator correspondence these degenerate composites represent the states that are 
holographically dual to the degenerate fermionic gas in the bulk theory. 

2.1.1 Degenerate fermionic operators 

Gravity in AdS^+i is holographically dual to a conformal field theory on S"^ -1 x time. The 
Hilbert spaces of the two theories are isomorphic, so any quantum state in the bulk is dual 
to a state in the CFT. The states of the CFT can be conveniently labeled by local operators 
on the plane via the state-operator correspondence 

lim §(z) I vac) = 1$) 

The energy of the state |$) is equal to A/ £ where i is the radius of S d ~ l and A the dimension 
of the operator $. The other quantum numbers of the state (like the angular momentum and 
R-charge) are also determined by those of $. 

In AdS / CFT one makes a distinction between single trace and multi trace operators. The 
first kind may be thought of as single particle states, while the latter are multi-particle states. 



8 However, as we will discuss in section [6j the N = 4 SYM may not be the best setup to realize the 
degenerate star. 
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However, this distinction is only well defined at large N, or more generally for certain classes 
of CFTs with a large central charge c. The reason being that 1/N effects will lead to a mixing 
between single trace and multi-trace operators. 

For certain strongly coupled supersymmetric CFTs the bulk theory is given by an AdS 
supergravity theory, which for large central charge can be treated (semi-) classically. The bulk 
theory contains bosons and fermions, and accordingly the boundary single trace operators 
may be distinguished into bosonic and fermionic operators. Whether a field is fermionic or 
bosonic is determined according to the spin-statistics theorem through its representations of 
the SO(d,2) conformal symmetry group. In particular, primary operators that transform 
in the spinor representation correspond to fermions, and have correlation functions that are 
anti-symmetric under the exchange of two identical operators. 

The fermionic nature of primary operators can in specific cases also be understood from 
their construction in terms of the basic fields of the underlying field theory. To give a con- 
crete example, let us consider the case of Af = 4 SYM. The spectrum of operators contains 
chiral primary fields that are holographically dual to the Kaluza-Klein modes of the type IIB 
supergravity theory on AdSsxS 5 . These chiral primaries form representations of the 5*0(6) 
R-symmetry group, and are arranged in short supermultiplets. The lowest component is given 
by the bosonic single trace operators 

Qhh...h = [ tr ( V 2 /3 • • • ^ symmctric 

traceless 

By applying a supersymmetry transformation one learns that the same supermultiplet con- 
tains fermionic single trace operators of the form 

K^ Jk = IQa*, * hh - h ] = [I^tr (A?0V" • • • <P h )} —c (2.1) 

traceless 

Here is an SO (6) gamma matrix, with spinor indices A, B and vector index I. We 
note that this operator is not primary with respect to the super conformal group, but it is a 
conformal primary operator with respect to its bosonic subgroup. The BPS bound implies 
that the field \& has conformal dimension k + \ for all values of the 't Hooft coupling A. 
The fermionic field \& also carries an SO (4) spinor index a, and a multi-index made from a 
symmetric traceless combination of vector indices. In the following we clearly wish to avoid 
using all of these indices. Indeed, except for the fact that these fields are fermionic, we have 
no need for any of these quantum numbers. This is because we will work either in a regime 
in which we ignore all interactions, or in an approximation in which we do the same except 
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for gravitjo All this will not be very relevant for what follows, but merely was meant to give 
a concrete example of a fermionic primary field. 

Let us now start with a fermionic single trace operator \I/ of dimension Ao- For simplicity 
we ignore the spinor indices of \I/ since this will not affect any of our following results qualita- 
tively. We emphasize that \& is not one of the fundamental field of the Lagrangian but rather 
a gauge invariant operator which has finite conformal dimension at strong coupling. We work 
in normalization where the 2-point function is 

(tt(x)tf (y)) = -. 

The 't Hooft large c factorization implies that for single trace operators the connected corre- 
lators are suppressed as 



which in turn implies that the field behaves effectively like a "free" field at large c. We 
can construct multi trace operators by taking products of many of these fields without having 
to worry about operator mixing with operators with a different number of traces. In a sense, 
at large c the fermionic operators generate a Fock space, just like free fieldio. This fact will 
be true as long as the number of operators, or more precisely the total conformal dimension, 
is small compared to the central charge c of the CFT. This point will become important when 
we discuss the possible role of interactions. 

We now turn to the multitrace operators that can be constructed from the operator 
$ '. From the fact that correlation functions are antisymmetric under the exchange of two 
fermionic primaries 

(... V(xMy)... > = -(... ¥(y)*(z)...) 
one derives that the operator product of two of these primaries is also antisymmetric 

V(x)iff(y) = -^(y)^(x). 



9 And even if we would try to keep track of R-charges, we can still construct neutral composite states 
by including all fermionic operators in a complete multiplet. The only complication this will lead to is 
a degeneracy factor equal to the dimension g of the SO (6) representation corresponding to the traceless 
symmetric s-tensors. 

10 Though it is not a standard free field in d dimensions, but rather a "generalized free field" [23]: it does not 
obey any linear wave equation since its conformal dimension is non-canonical but nevertheless its correlators 
factorize. 
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Normally one defines the normal ordered product : ty 2 (x) : by taking the limit x — > y after 
subtracting possible divergent terms. For fermionic fields the result is zero, due to the anti- 
symmetry. This is of course just a reflection of the Pauli exclusion principle. Therefore, to 
construct a nontrivial operator out of two fermionic primaries we have to include a derivative. 
This gives operators of the kind : \I/<9j\I/ :, where di denotes the derivative w.r.t. X{. This 
operator is the operator of lowest conformal dimension that can be found in the regular part 
of the operator product expansion of ^(x)^(y) [3 For a <i-dimensional CFT there are d types 
of these operators. 

We want to continue with the construction of multi-trace operators of the lowest possible 
conformal dimension made out of three, four, all the way to a very large number of copies 
of the basic fermionic operator Each time we have to consider the fact that the next 
operator needs to be anti-symmetrized with the previous ones. So we have to act again with 
derivatives every time making sure that the combination of derivatives has not been used 
before and also that we are using the smallest number of derivatives that are necessary. If 
we have constructed an operator of this type with Np insertions of the basic operator ty, we 
can inductively define the degenerate operator with Np + 1 insertions as the operator with 
lowest conformal dimension appearing in the regular part of the OPE between the composite 
with Np insertions of \l/'s and one additional This definition is unambiguous at infinite c. 
Continuing like this one notices that one is filling up " shells" , very much like an atom, where 
each shell is labeled with the number n of derivatives that are being used. 

In order to end up with a rotationally symmetric composite operator it is natural to 
consider the "noble elements" with a completely filled last shell. We will denote the number 
of derivatives in the last shell by np, since it plays the role of a Fermi level. In this way we 
arrive at the following form of the degenerate composite operators 

*=*Y[a i *Y[d i d j v n didjdkV n d h d i2 ...d in y (2.2) 

« {m} {i,j,k} {h, i 3 , ...i„} 

where one takes in each shell the product over all possible n-tuples made from d derivatives. 
We note that the composite operator only depends on the choice of the number np. It does 
not carry any other indices, and hence the corresponding state is unique. In other words, we 
do not have to deal with counting states: the total entropy is equal to zero. It should also be 
obvious why these composite operators are called "degenerate". 

In the bulk interpretation each insertion of the operator \1/ represents the addition of a 



This statement is precise in the c — > 00 limit. 
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particle. Some simple combinatorics show that a completely filled shell made from fields with 
n derivatives contains 

fn+d-1 
\ d-1 

different fields, and hence corresponds to that many particles in the bulk. The total number 
of fields (or particles) is found by summing over all shells 

^ = £ ("d-i 1 ) = (" F d d ) ^ 

Next let us now determine the conformal dimension of a degenerate composite operator $ 
with np filled shells. The conformal dimension of the operators in the n-th shell is n+ A , since 
each derivative adds one to the conformal dimension. Hence the total conformal dimension is 

A = g(„ + A„) = - (Xl) + *° ^T) ^ 

The second term is proportional the the total number of fermions and represents the contri- 
bution of their "rest mass" to the energy, while the first term is the "kinetic energy" coming 
from the derivatives. 

In this discussion we have ignored so far the spinor index of the field. This can be treated 
in a similar way as the electron spin in atomic physics. It simply tells us that each state 
can be occupied by a number of fermions, namely as many as the dimension of the spinor 
representation. So if we would have included it, we would have to write an additional product 
over a which would have multiplied the counting formulas for Np and A by the dimension 
of the spinor representation. Similarly, if we would have considered a number of species g of 
fermions, these formulas have an overall factor g. 

Apart from these harmless simplifications we would like to emphasize that our discussion 
so far has been exact for all values of the coupling of the gauge theory, as long as c — > oo. We 
have assumed that the basic building block \1/ of our operators has conformal dimension A . 
At the moment we do not care whether this Ao is the same as the conformal dimension of \I/ 
at weak coupling. For us A can be considered as an input parameter which is determined 
by the (possibly strongly coupled) dynamics. But even this issue could be settled by taking 
$ to be a chiral primary in a supersymmetric theory. On the other hand the operator <3? 
defined in (12.21) is definitely not supersymmetric, even if \1/ is, since the derivatives increase 
the conformal dimension without adding R-charge. The reason that the expression (12.41) is 
still reliable is based on the 't Hooft factorization at large c: the conformal dimensions of 



11 



multi-trace operators are equal to the sum of the individual conformal dimensions up to - 
corrections. This statement is based on the planar expansion and does not depend on the 
value of the 't Hooft coupling. Hence in the limit c — > oo our formula (j2.4p is exact. The 't 
Hooft factorization also guarantees that at large c we do not have mixing of the operator $ 
with other operators with the same quantum numbers. These statements are reliable if we 
take c — > oo keeping other quantities (like A and nf) fixed. 

Finally, let us finish this section with a comment on the SO(d) decomposition of the 
states in the n-th shell. These states are given by symmetric n-tensors, but since they are 
not traceless they do not form an irreducible representation. The dimension of the irreducible 
representation of the symmetric traceless n-tensors is given by the number of symmetric n- 
tensors minus the number of symmetric (n—2) -tensors. In this way we find that the states in 
the n-th shell decompose into a sum of representations corresponding to traceless symmetric 
(n — 2k) -tensors, with k an integer that runs fronO to [n/2] . Just as a check one verifies 
that 

n + d — 1 
d-1 



[«/2] , 

' n-2k + d-l\ f n-2k+d-3 



£ 

k=0 



d-1 / V d-1 



(2.5) 



2.1.2 "Free field" mode expansion 



The state-operator correspondence is a convenient way to describe states in a CFT. The 
relationship with the corresponding bulk state becomes more manifest, however, in a slightly 
different representation of the states. Instead of creating the state at r = — oo one can 
represent the same state by making use of a mode expansion of the primary field on the 
cylinder S^ -1 x R. So we go back to Minkowski signature. At large c the primary operator 
can be expanded in modes [24J on S as 

*(t,n) = ^a nAm yK^-i)e^ + <, im y;(^-i)e-^ (2.6) 

n,l,m 

where the energies satisfy the condition E n \ = (A + n)/£, n = 0,1,2..., determined by 
the representations of the conformal group. The expansion (12.61) has to be understood as a 
definition of the operators at n) i )m , /3 n ,i,m- Using the 2-point function of the field \I/ and the 
fact that its correlators factorize at large c we find that the operators cx nt i >m , p n ,i,m satisfy the 
algebra of fermionic ladder operators 

{® n ,l,m, <^ n i i> m >} ~ <W<V^m,m' > {Pn,l,m, 0tf y m '} ~ ^nn'^U'^m^m' (2.7) 



12 Here with [ .. ] we denote the integral part of a possibly fractional number. 
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and all other anticommutators are zero. The vacuum state satisfies the condition a nj i jm \vac) = 
/3 n ,i,m\vac) = 0. So we can use these modes to build a multi particle state from the vacuum. 
In this way one finds the state |$) corresponding to the operator (12.21) has an alternative 
representation as 

i*> = n a iiJ vac ) ( 2 - 8 ) 

n,l,m 

However it should be stressed that the expansion (12.61) and the commutation relations (12.71) 
are only true in the c — > oo limit, keeping other quantities fixed. If the operator \I/ is 
inserted between states whose conformal dimension is of the order of the central charge, then 
factorization may not hold any more and the commutation relations (12.71) will be modified. 
Hence the repeated use of the modes of the operator \1/ is only allowed in the limit of large 
central charge where factorization is reliable. 

2.1.3 Density of states 

Let us now introduce the total mass M, the fermion mass rrif, and the Fermi energy ep via 

A ° m A n F + A 

where i is the radius of the boundary S d ~ x . We express the number of particles Np given in 
(12. 3p and the mass M in terms of the Fermi energy ep and the fermion mass m/ in the limit 
that np = (ep — mf)t and Ao = rrifl become large. The number of fermionic particles may 
be expressed as 

Np = j l {e F -m f ) d (2.10) 

where we ignored all terms that are sub leading in rip. For the total mass we use the expression 
(12 ,4p for the total conformal dimension A. Keeping the leading contributions of both terms 
we find 

M = ^(e F - m f r + ^(e F ~ m f ) d (2.11) 

It is important that we keep both ep as well as mj in these expressions, since we want to 
keep the ratio of ep and rrif finite. The expressions (I2.10p and 12. 11]) only hold in the specified 
regime, since we have not taken into account any interactions, gravitational or otherwise. 

We will be interested in the limit of large np (and large A ), so it is more convenient to 
rewrite the sums in the expressions f ]2. 31) and ( 12. 4p as integrals. This gives 

N{e F )= / deg(e), M{e F ) = / deeg(e), (2.12) 

J mj J rrif 
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where the density of states g(e) is given by 

This expression is valid in the strict large c limit. 

Notice that this density of states differs qualitatively from that of a free fermion in d 
spacetime dimensions. For the latter the density of states grows like e d ~ 2 . In other words, 
while the fermionic operator \1/ is defined in a d dimensional CFT we have found that the 
density of states created by it grows like that of a d + 1 dimensional gas. Of course this is 
related to the fact that the operator \1/ is not really a fundamental free field obeying a linear 
wave equation, but rather a "generalized free field" which naturally lives in one additional 
dimension. In a sense this dimension is the scale in the conformal field theory. 

2.2 Quantum states in the bulk 

The fermionic single trace operator ^ corresponds to a fermionic field ip in AdS spacetime. 
On-shell this field satisfies a field equation which in principle depends on the spin of the 
field. For spin one-half it would be the Dirac equation, or for higher spin one would get the 
Rarita-Schwinger equations and its generalizations. But since eventually we are interested 
in taking the limit of a large number of particles and making contact with a hydrodynamic 
description, the distinction between the different type of fermionic particles and field equations 
becomes irrelevant. Therefore, to simplify our discussion we will use the partial waves of the 
Klein-Gordon equation, which are easier to deal with and closer to the standard knowledge^! 
The global AdS^+i metric is given by 



ds 2 = -[l + -j dt 2 + dr 2 + r 2 dn 2 d _ v (2.14) 

where £ is the AdS radius. 

We will assume that the field ip satisfies the rf+l-dimensional Klein Gordon equation in 
AdS space, which in full detail looks like 



1 + i I d 2 + r l - d d r ((l+ r -p)r d - x d r ) + r- 2 A n + m) 



ip(t,r,n) = 0, (2.15) 



13 The fact that the Klein-Gordon equation usually describes bosonic particles does not pose a problem, 
because we can implement Fermi statistics by hand and take the degeneracy of the spin degrees of freedom 
into account afterwards. Again, this is very analogous to the standard treatment of the electron spin in atomic 
physics and is a reliable approximation in the limit of large number of particles. 
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where A^ denotes the laplacian on S 1 ^ 1 . The mass of the field ip and the conformal 
dimension of the operator A are related by the formula A = & + ^/f + ( mf l)*. For 
large values of my we can approximate A ~ rrifi, which of course was the reason for the 
identification in (12. 9p . Since we are interested in states which can be considered as finite energy 
excitations above the standard vacuum of the theory, we only focus on the normalizable modes 
in AdS. We then find that the partial wave solutions of the Klein Gordon equation take the 
general form 

Vv,m(*, r, n) = e -*(™+ A °)^ / n _ 2M ( r ) Y l>m (Q). (2.16) 

Here / n ,«( r ) * s some hypergeometric function [25], whose precise form will not be important 
for us. The key point is that both the bulk and the boundary fields have the same quantum 
numbers. This makes the correspondence straightforward. To be precise, the relationship 
with the operator \1/ in the CFT may be expressed by 

#(t, n)\vac) = lim e A ° r/e tp(t, r, Sl)\vac) (2.17) 

r— >oo 

This implies that the creation modes that we defined in the CFT are indeed directly related 
to those constructed out of the bulk fields. The Hilbert spaces are therefore completely 
identified, at least in this free field limit. This also implies that the composite operator that 
we constructed indeed may be identified with a multiparticle state, whose complete wave 
function may be represented in the bulk as a Slater determinant in terms of the modes given 
in (12TT6D . 

i*) = n a UJ vac ) ( 2 - 18 ) 

n,l,m 

2.3 Hydrodynamic description 

In this section we will give a hydrodynamic description of the multi-particle state as a free 
degenerate Fermi gas in AdS space. The Fermi gas is found to be confined to a spherical 
region with a radius determined by the ratio of the Fermi momentum and the fermion mass. 
This fact is due to the gravitational potential well of the AdS space and not due to self-gravity. 

2.3.1 The equation of state 

At first one may think that the equation of state for a degenerate Fermi gas in AdS space is 
different from that in flat spacetime. For instance, the spacing of the lowest energy levels are 
determined by the AdS radius, and differ from that of flat space. But when the number of 
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particles within one AdS volume is very large, one can safely ignore the effect of the curvature 
radius, and use a local description of the gas as if it lived in flat spacetime. As we will see this 
approximation reproduces the exact results of the previous subsection in the limit of large 
number of particles. 

To make this point more clear, consider a degenerate free Fermi gas in a box of size L in 
d spatial dimensions. The Fermi momentum kp is related to the chemical potential \x via 



H = yJk F + m 2 f . (2.19) 



At zero temperature all states satisfying the bound 

2ii\n 



L 



< kp (2.20) 



are being occupied. Semi-classically we can count states by calculating the volume of phase 
space. Hence, the particle number density n is equal to the volume of the d-dimensional 
Fermi sphere: 



n= ^zlZJL. (2.21) 

Here V^-i is the area of a unit (d— l)-sphere. 

These are the usual relations that hold in flat space. But as long as the number of particles 
per AdS volume is very large, the Fermi momentum will be much larger than the inverse AdS 
radius t. Therefore, in order to determine the equation of state one can take the size L of the 
box very small compared to £. This means that the flat space relations given above also apply 
locally in the curved AdS background. The only difference with a translationally invariant 
flat space situation is that in this case one should expect that the Fermi momentum kp and 
chemical potential \i depend on the position in the AdS space. 

The energy density and pressure are given by the standard expressions 



V, 



(2 



7T 



"F 



d_1 / At. Ud-l U2 i ™2 







dk k d - L Jk 2 + m 2 f , 



v = V ^~ 1 f k F dk kd+ — - (2 22) 

and obey the standard thermodynamic relations 

dp — pidn, p + p = fin. (2.23) 
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One small comment about the energy density p will be useful. Namely, it can be written as 

dn 

de 



i ep dn , 

p= I e — de (2.24) 



Intuitively, the validity of this identity should be obvious: one simply adds all energies e with 
a weight given by the local density of particles with that energy. Similarly the pressure can 
be written as 



p = I n(e)de 



2.3.2 Energy momentum conservation. 

In the hydrodynamic limit we expect to be able to describe the free Fermi gas purely in terms 
of the local energy density p, the pressure p, and a local velocity field w M normalized so that 
u^u^ = — 1. The energy momentum tensor is given in terms of these quantities as 

Tftu = (p + p)tyu„ + pg^ (2.25) 

We will first consider the non-rotating degenerate state and postpone the discussion of the 
rotating state to section 12.41 Without rotation the energy density and pressure will only 
depend on the radius r. In this section we are still ignoring the effect of self-gravity of the 
gas. However, we still have to take into account the effect of the curved AdS background. This 
is done by imposing conservation of the energy momentum tensor. For a static configuration 
it is the equivalent of imposing hydrostatic equilibrium. Anticipating possible generalizations 
let us write the metric in the general form 

ds 2 = -A{rfdt 2 + B{rfdr 2 + r 2 dn 2 d _ v (2.26) 

The Fermi gas is static, and hence we take ut = A(r) and all other components of u equal to 
zero. The only non-trivial equation is V M T Mr = 0, which becomes 

| + + = 0. (2.27) 

This equation is surprisingly easy to solve. By making use of the identities (I2.23P one easily 
verifies that (12.271) is satisfied when the chemical potential obeys 

Mr) = ^r y (2.28) 
17 



where at this stage ep is an arbitrary constant. We conclude that the radial dependence of 
the chemical potential is simply due to the gravitational redshift. With some hindsight this 
result could have been anticipated: it is an analogue of the familiar Tolman relation for the 
temperature, which is known to hold for any static gravitational field. 

So far, this discussion has been general in the sense that we have not used the explicit 
form of the AdS metric. But let us now go back to the specific case of the AdS space time, 
and investigate its energy and pressure profile in more detail. For the AdS metric ( I2.14p the 
relation (I2.28P becomes 

"W = 7=W' (2 - 29) 

Now that we have determined the radial profile p(r) of the chemical potential we in principle 
also know the value of the energy density p(r) and the pressure p(r) through the relations 
(I2.22p . It is clear from this result that the chemical potential decreases as one goes out to 
the boundary. It can not decrease arbitrarily, however, since the energy per particle is always 
larger than the mass rrif. Namely, when // = rrif the Fermi momentum kp vanishes and so 
do the energy density and pressure. We conclude therefore that p(r) and p(r) will go to zero 
at a finite value R of the radius, namely when 

ljl{R) = to/. (2.30) 

It will be interesting to re-express the value of the radius R in terms of the Fermi momentum 
kp(0). Combining the equations (I2.29f) and H 2 . 3 1) together with the standard relation between 
kp and ep gives 

R= k -^L (2.31) 

TO/ 

Hence, the size of the star in AdS units is a direct measure for how relativistic the fastest 
particles are in its center. When the size of the star is of the order of the AdS radius i these 
particles are mildly relativistic, while they have ultra-relativistic velocities when the star is 
much larger than I. It is important to note that the Fermi gas has a sharp edge only when 
the fermions are massive. For massless fermions the equation of state is simply that of a 
relativistic gas, and the star will spread out with a dilute "tail" towards the boundary of 
Ad£& 
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Of course when the density gets too low the hydrodynamic approximation breaks down. 
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2.3.3 Total mass and particle number. 

The degenerate star-like object that we described in the previous subsections is the holographic 
dual of the boundary state |<&) constructed in section 2.1. As a non-trivial check on our 
hydrodynamic description one can verify that the particle number and mass agree with the 
boundary calculations. We observe that the constant ep equals the chemical potential defined 
with respect to the time t. This means that it can be identified with the Fermi energy ej? of 
the CFT 

ef T = e AdS (2.32) 
The total particle number in the bulk is given by the integral 

r R 

Np dS = Vd-i / drr d ~ l B(r)n{r). (2.33) 
Jo 

The factor B(r) comes from the volume form on the spatial section. Similarly, one obtains 
the mass M AdS by integrating the energy density. Because M AdS is defined with respect to 
the time coordinate t one has to include a redshift factor A(r). For the AdS metric ( I2.14p the 
redshift factor A(r) cancels the measure factor B(r), and hence 

r R 

M AdS =Vd_! drr d - l p(r). (2.34) 
Jo 

To calculate Np dS and M AdS one has to insert the equations (I2.2ip and f!2.22p for the particle 
and energy density and include the redshift effect (12. 28ft into the chemical potential. The 
resulting integrals can be performed analytically and (perhaps not surprisingly) precisely 
reproduce the results (I2.12p and (I2.13P in the limit of large number of particles. 

2.3.4 A comment on the validity of the hydrodynamic description 

In general, a basic assumption of the hydrodynamic approximation is that the continuous 
system under investigation is in local thermodynamic equilibrium. More precisely the hy- 
drodynamic description is an expansion in (space-time) derivatives of local thermodynamic 
quantities. As such, hydrodynamics is valid when the typical variation length is much larger 
than the microscopic scales of the fluid. In particular it should be much larger than the mean 
free path. In our case the "fluid" is actually a free gas, so the mean free path is infinite or 
at most bounded by the AdS scale. This means that strictly speaking we are not allowed to 
use hydrodynamics. In particular if we had considered time dependent configurations, then 
there would clearly be a problem since our fermions are non-interacting and there would be 
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no restoring agent to bring them back to equilibrium^ . However for our static, spherically 
symmetric situation the results derived by the hydrodynamic method are reliable. This is 
indicated by the agreement of (I2.10p . (12. lip with (I2.33p . (I2.34p and could perhaps be made 
more quantitative along the lines of [26H28] . 

2.3.5 Radial Reconstruction 

In the limit of a large number of fermions the hydrodynamic description in the bulk is more 
convenient than keeping track of the exact quantum state (I2.18p . We saw in sections 12.11 and 
12.21 that at the level of multiparticle quantum states we have a one-to-one correspondence 
between the boundary and the bulk. This means that we should be able to construct the 
analogue of the "hydrodynamic description" in the boundary CFT. In other words, we would 
like to understand what is the CFT meaning of a hydrodynamic description of the multitrace 
operator ( 12. 2p . 

The hydrodynamic description of the gas in the bulk is partly encoded in the radial density 
profile p(r) given by (12.221) and f!2.28[) . How is p(r) related to the operator (12. 2p ? The radial 
direction r does not have a direct meaning in the CFT, though in a sense it corresponds to 
the "scale" of the excitation. One way to probe the scale direction is to consider 2-point 
functions of operators in the CFT, evaluated on the state (12. 2p . The separation of the two 
insertions of operators on the boundary controls how deeply in the bulk we are probing: when 
the two operators are close to each other we are probing the region near the boundary, while 
at larger separations of the boundary operators we probe the interior and in particular the 
region where the fermionic gas lives. Hence if we compute the 2-point function 

in the limit where the number of particles in is large, we should be able to relate the 
radial profile p(r) to the (x — y) dependence of this correlator. 

Computing this correlator should be straightforward in the free c — > oo limit and it would 
be interesting to perform this computation, as this would be the first step before considering 
the same problem including backreaction. We hope to revisit this question in future work. 



15 Once we consider 1/N corrections the particles will become interacting and they will of course reach 
equilibrium. 
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2.3.6 Fermions at finite temperature 

It is not difficult to generalize the story to the case where the fermions are at non-zero 
temperature. Let us assume that the fermions are in an ensemble characterized by the average 
occupation number £(e) for the one-particle states which are multiplied together to make the 
analogue of the multi-trace operator (I2.2[J ^. In the zero temperature case the occupation 
numbers are given by the zero temperature Fermi-Dirac distribution £(e) = Q(e F — e), where 
G is the unit-step function. This distribution corresponds to the radial profile ( 12. 291) of the 
local chemical potential in the bulk. We denote the corresponding density and pressure profiles 

by Pe F (r),Pe F ( r )- 

What happens if we consider a more general distribution £(e) for the average occupation 
level of the single-trace operators ? Notice that formally we can write any function £(e) 
defined for e > as a linear combination of O-functions in the following way 

poo 

£(e) = -/ de F e(e F )Q{e F -e) (2.35) 
Jo 

which can be proven by a partial integration, assuming that £ vanishes at ep = oo. Since the 
fermions are non-interacting this implies that the bulk profile can be written as 

/•OO /*O0 

Pd r ) = ~ de F ^'{e F )p eF {r) , p € (r) = - / de F £'(e F )p tF (r) (2.36) 
Jo Jo 

In particular for a thermal distribution in the grand-canonical ensemble at temperature j3 and 
chemical potential \x we have 

£( e ) = —— (2.37) 

Using (I2.29p . f |2.36|) and (I2.37j) we can compute the radial profile of the density (and pressure) 
corresponding to a finite temperature fermion gas in AdS, before we consider any backreaction. 
Let us notice that while the equations (I2.36P seem to work for any choice of £, the profile 
in the bulk will have a meaningful description in terms of a density profile p%(r) only if £ is 
sufficiently smooth. 

It would be interesting to compute the entropy of a thermal ensemble from the fluid 
configuration in the bulk. For large quantum numbers this entropy should agree with the 
exact boundary computation of the entropy in terms of the quantum states of the single-trace 
operators. 



16 For simplicity we assume that we are in an ensemble where the average occupation level £(e) depends only 
on the energy of the state. This assumption guarantees that the configuration will be spherically symmetric 
in the bulk. 
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Finally, from the relation ( I2.36P it should be clear that the mapping between perturba- 
tions of the bulk density profile and the boundary is non-local, as expected from the general 
framework of AdS/CFT. 



2.4 Rotating Star 

As a straightforward generalization and as a further test of the hydrodynamic approximation 
we will now try to incorporate rotation into the system. For simplicity we focus on AdSs and 
we take the metric to be 

ds 2 = - + dt 2 + (l + dr 2 + r 2 {d9 2 + sin 2 6d<f) 2 + cos 2 Od^ 2 ). (2.38) 

The range of the variable 9 is [0, 7r/2]. In order to put angular momentum into the system 
we need to take a rotating fluid, and there is an obvious question whether or not the rotation 
is uniform. In four dimensional flat space it is claimed in [29], [30] that uniform rotation 
minimizes the mass of the star. We are going to try to do the same thing for the five- 
dimensional AdS case. 

The energy momentum tensor still equals 

Tpu = (p + p)Ufj,u u + pg^ (2.39) 
where now p and p will be functions of both r and 9, and 

= u-\i,o,o,n^,n^ (2.40) 

where ^ 

u = ( + i) ~ r2 sin2 * n * " r2 cos2 1 2 ' (2 ' 41) 

Conservation of T implies 

= U 2 d r p + £(p + P ){l-en% cos 2 9 -£ 2 0|sin 2 9) 

= U 2 d e p + r 2 (p + p)(nl-nl) sin 9 cos 9. (2.42) 

These equations are very simple, they are equivalent to 

= d r p+^-(p + p) 

= d e p+^(p + p). (2.43) 
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We can now immediately solve for the star, the solution simply reads 

,M = ^ (2.44) 

If Q^, depend non-trivially on r, 9 this simple solution is no longer correct. 

Given the exact solution, we wish to write expressions for the mass and two angular 
momenta. These are all related to Killing vectors ^ and the conserved charges boil down to 

Q [£] = 4tt 2 J f " T m0 ( X + ^2 ) rS sin 6 cos 9 dr d0 - ( 2 - 45 ) 
In particular, we obtain for the mass 

M = Att 2 J r 3 sin 9 cos 9drd9 (-p + 1 (p + (2.46) 
and for the angular momenta 

J^ = -4vr 2 y r 3 sin cos flrfr rfg ^!^_ % + p ). (2.47) 

and 

= -An 2 j r 3 sm9cos9drd9^-^^-(p + p). (2.48) 
These are rather difficult expressions, which simplify quite a bit if we put 

= ^ = |. (2.49) 

One can then explicitly do the relevant integrals, and we obtain 

M _ & (fg ~ "V) 4 (4e F + m f - 5m f u 2 ) 

M ~ i2o (T^r { } 

and 

J " 60(1-^)3 • V- 51 > 
Let us now describe the rotating degenerate state in the boundary CFT. For this purpose 
we choose a U(l) subgroup in one of the two SU(2) factors. To be specific, let us consider 
a rotation associated with, say, a simultaneous right handed rotation in <pi and <p 2 - The 
corresponding quantum number is 

j = mi + m 2- (2.52) 

Single particle states are characterized by their angular momentum I and energy E = A + n. 
When n is fixed but / is left free j takes the following values 

-n<j<n, (2.53) 
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but remains even (odd) when n is even (odd). 

We describe the rotating degenerate state by introducing an angular velocity u which acts 
as a "chemical potential" for the quantum number j. This means that for a given Fermi 
energy e F all states are being occupied which obey the inequality 

n + A - uj < e F I. (2.54) 

Thus the effect of the rotation is to tilt the Fermi surface to favor positive values of j. 

To count the total number of particles contained in this rotating degenerate many particle 
state we have to add all one particle states labeled by n and j that satisfy the inequality 
(I2.54p . As a preparation, let us count how many single particle states exist for given values of 
n and j. This is most easily done by making use of a generating function for the symmetrized 
tensor products of the vector representation of SO (4). One can easily convince oneself that 
this generating function is given by 

1 (2.55) 



(1 -qz) 2 (l -qz- 1 ) 2 ' 

Here q keeps track of n, while z keeps track of the j quantum number. Picking the term 
proportional to q n z^ gives the following degeneracy 

A/-(n,j) = ^((n + 2) 2 -j 2 ) (2.56) 

Note that for even (or odd) n the highest number of states occurs at j = (or j = ±1). For 
j = ±n one finds precisely one state, as expected. 

This result can now be used to count the number of particles in the rotating degenerate 
state for a given chemical potential e F . We will do so in the limit of many particles, so that we 
can replace the summation over the quantum numbers n and j by integrals. In fact, for this 
purpose it is convenient to replace n + 2 simply by n. This leads to the following expression 

N(e F ,cj) = l j Af(n,j)dndj. (2.57) 
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where T> tP}UJ denotes the integration domain bounded by (I2.54p . (I2.53|) and n > 0. The 
factor of 1/2 is due to the fact that the values of j are only even or odd. The integrals are 
straightforward but somewhat tedious. After some work we find 
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We see that for a given chemical potential the number of particles increases due to the 
rotation. Note further that in the limit u — > ±1 the number of particles diverges. This can 
be understood from the fact that in this case all states with j = ±n are being counted for 
arbitrary values of n. 

The mass and angular momentum can be obtained in a similar way as integrals over n 
and j over the domain T> ep U] . For the mass one has to sum up the energy e = (n + A ) /£ for 
all particles contained in the rotating degenerate state. This leads to the integral expression 



M = hl (ra + Ao) A/ ' (n ' j) dn dj - 



(2.59) 



Performing this integral is again straightforward. The final result can be organized in the 
following form 

^ Uf ~ m o) 5 rani 4 Up — mo) 4 

This result should be compared to the case without rotation given in (12.111) . The last term 
represents the contribution of the mass m for all the particles, whose number is given in 
(I2.58p . while the first term gives the sum of all the n quantum numbers. We note that the 
modification in the mass due to the rotation is again given in terms of inverse powers of 
(1— u 2 ), and diverges for u — > ±1. 

Finally, let us come to the angular momentum J. The total angular momentum carried by 
the particles in the rotating degenerate state is represented by again a tedious but fortunately 
straightforward integral 

J= \t\ ^(nj)dndj. (2.61) 

It yields the following result 

J= *<■*'- m ? u (2.62) 
60(l-w 2 ) 3 V ; 

Note that \J\ < M in general, and only when u — > ±1 the ratio \J\/M — > ±1. 

The results (I2.60p . (I2.62p for the mass and angular momentum computed in the conformal 

field theory agree with the hydrodynamic computations (12.501) . (I2.5ip . 

3 INCLUDING INTERACTIONS 
3.1 The scaling limit 

The discussions of the previous section are accurate if we take the c — > oo limit keeping A, Np 
and A fixed. In that limit the system is free as interactions can be ignored due to 't Hooft 
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factorization. However what we have constructed so far is a trapped fermionic gas rather than 
a star. To turn the gas into a star we have to include gravitational backreaction. This can be 
achieved by scaling Ao and Np (and as a result A) appropriately at the same time that we 
send c to infinity. 

To determine what is the precise scaling limit we should take, we consider the picture in 
the bulk. The first condition we impose is that the radius of the star is of the order of the 
AdS scale £ Assuming that equation (12.31 1 ) is a good estimate for the order of magnitude 



of the size of the star even after backreactiorj 18 !. we see that this condition implies 

hlL = ^L = fixed (3.1) 
mf A 

The second condition is that we want to have appreciable gravitational backreaction. This 
can be translated to the statement GM£ 2 ~ d ~ 0{1). Newton's constant is related to the 
central charge by G ~ c~ x l 1 . So we find that the conformal dimension of the composite 
operator must satisfy 

— = fixed (3.2) 

c 

where we used the identification A = Mi. From equation (12. 4p and (13. ip we find that for this 
to be true we have to take 

n F , A ~ c~ 
In this limit the number of particles ( I2.3P goes like 

d 

Np ~ c d +! 

In other words we want to study a specific class of composite multitrace operators in the limit 
c — > oo , A = ec where e is kept fixed. When £ < 1 the system is weakly interacting. We 
want to turn on e up to values of order 1 where interactions become importantly. 

Since in this limit the number of particles becomes large, it is conceivable that the 't 
Hooft factorization might fail in a drastic way. In terms of Feynman diagrams in a gauge 
theory this can happen because we are scaling the number of external lines of a diagram 
with N so the usual iV-counting rules may not apply. Let us describe the problem in some 
more detail: consider the correlation functions of multitrace operators whose size we keep 



17 By this we mean that as c — > oo the radius of the star scales like as R ~ c°£. 
18 This assumption will be proven to be self-consistent after we actually compute the backreaction. 
19 Let us note that even though in our limit the particles become very heavy (since their rest mass scales 
like £ _1 c^+ T ), they are parametrically lighter than the Planck mass which scales like ^ _1 c 3zrT . 
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fixed in the large N limit. As far as interactions between distinct single trace constituents 
are concerned, the dominant diagrams are the disconnected ones, which correspond to the 
free limit discussed in section |2j Connected planar diagrams involving interactions between 
various single trace components of the multitrace operator are suppressed by powers of 
relative to the disconnected ones and non-planar diagrams by even higher powers of j*. Now 
let us imagine scaling the size of the multitrace operators as we send N — > oo. In order to have 
nontrivial interactions at large N we want to scale the size of the multitrace operators in such 
a way as to compensate for the -k suppression of connected planar diagrams (connecting the 
various single trace components of the multi-trace operators). Our goal is to work in a limit 
where these leading connected planar diagrams are of the same order as the disconnected ones. 
The danger now is that since in this limit the planar connected diagrams became as important 
as the disconnected ones, then the same could happen for non-planar diagrams. This would 
imply that effects of quantum-gravity would become important in the bulk. However we will 
argue below that this is not the case and the planar diagrams remain the important ones. 

It is perhaps more intuitive to address this question from the bulk point of view. In theories 
with a gravitational dual the planar interactions on the boundary are mapped to tree level 
supergravity diagrams in the bulk. The c — > oo limit corresponds to sending G — > in AdS 
units. In this limit the tree level interaction between any two particles (whose mass is kept 
fixed) will vanish. In our double scaling limit we are taking the number and masses of particles 
to be large so that even though G — > the total backreaction remains finite. Moreover the 
various factors contribute in such a way that the backreaction receives contributions from 
tree level diagrams only, as we explain in the next subsection. This implies that the planar 
approximation remains reliable in this limit. 

3.2 Interacting Quantum states in AdS 

Let us now consider the many-body state (12. 18[) in the scaling limit of the previous subsection, 
in which we have to take into account the gravitational interactions between the fermions. It 
is intuitively clear that in the limit of large number of fermions per unit volume, the most 
efficient way to analyze the interactions is by working in the hydrodynamic description of 
section 12.31 but this time incorporating the gravitational backreaction of the fermionic fluid. 
This approach will be pursued in subsection 13.41 For now we will consider the interactions 
from a more microscopic point of view, in order to understand the validity of approximations 
involved in going from the exact to the coarse grained hydrodynamic description. 
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A first simplification is that in our scaling limit (and for theories with a semi-classical 
gravity dual) we can ignore processes which change the number of fermions. In a relativistic 
system the number of particles does not have to be conserved. If we put a certain number 
of fermions in AdS and then turn on interactions it is not guaranteed that their number 
will remain constant. Unless the fermions are protected by a conservation law, there will 
be processes which can change their number, as well as processes which produce all kinds of 
other particles such as gravitons etc. In principle all of these have to be taken into account. 
However in our scaling limit the fermionic gas becomes effectively long lived i.e. the rate 
of creation and annihilation of particles becomes negligible. This is due to the fact that at 
large c interactions are suppressed. In principle this suppression could be compensated by 
large combinatoric factors due to the large number of particles, but as we will discuss in more 
detail in section [BJ these factors work out in such a way that the processes which change the 
particle number become subleading at large c, relative to processes responsible for the total 
gravitational backreaction of the fermions. Because of this, in the rest of this section we will 
only focus on interactions which do not change the number of particles, i.e. on processes with 
only Np incoming and Np outgoing fermions and no other external particles. 

Second, we will assume that the only low lying fields in the bulk are the graviton, de- 
scribed by the Einstein-Hilbert action, minimally coupled to a fermionic field if) of mass rrif 
3 We will assume that there are no other types of particles or forces in the bulk. Clearly 
this is a very drastic approximation which however captures the qualitative behavior of the 
phenomena we want to study without unnecessary complications. It would be very interesting 
to study further whether a specific AdS/CFT setup can be found where such a simple bulk 
Lagrangian can be realized (at low energies), or at least a setup in which it correctly captures 
the qualitative physics of degenerate fermionic states. We present some first discussions along 
these lines in section [61 

We now want to see that in our scaling limit the - suppression of interactions combines 
with the enhancement from the large number of particles in such a way that only the tree 
level interactions (or equivalently planar interactions) become important. To illustrate this 
point we focus on one particular physical observable as an example: we will try to estimate 
the total energy of the fermions including their gravitational interactions. So let us consider 

20 There may be interaction terms for the field ip, but if their strength is controlled by c in the standard way 
(i.e. a vertex with n fermions is suppressed by a factor of c~) then the qualitative results of this section 
remain the same. For simplicity we will ignore such terms in what follows. 
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a bulk state |$) corresponding to Np fermions in the free theory. It satisfies 



H m = 




where H is the Hamiltonian for a free fermionic field in AdS^+i- When we include interactions 
we have to use the full Hamiltonian H = Hq + Hi nt where Hint includes the coupling between 
fermions and gravitons and the self-interactions of the gravitons and Hq has to be extended 
to include the quadratic part of the Einstein-Hilbert action for propagating gravitons. The 
corrected energy A/£ can be computed by an expansion in Feynman diagrams in the bulk. 

For example, consider two fermions in AdS with energies Ei, E 2 and wavefunctions ip\ (x) , ip2 (x) . 
To leading order the total energy is the sum of the energies. To compute the first correction to 
this energy we have to compute the tree level gravitational interaction between them. In the 
Born approximation this is given by the diagram shown in figure 1, where we have to use the 



Figure 1: Tree level gravitational interaction between two fermions. 



graviton and fermion propagators in AdS^+i- In our scaling limit this diagram is proportional 
to 

„mim2 _4=i 
G ~ c d + 1 

so it goes to zero. Higher order diagrams between two fermions are suppressed by further 
powers of -, so they all become negligible in the large c limit. Hence if we only had a finite 
number of particles the backreaction would be unimportant. 

If we have Np fermions then to compute the correction to the energy (or the conformal 
dimension of the dual operator) 5 A = A — A we have to sum over many kinds of diagrams 
as shown in figure 2 o, where as we explained before we only consider diagrams with the 
same number of fermions in the past and future infinity. Due to the non-linearities of general 
relativity we not only have 2-particle interactions but also n-particle ones for all n as depicted 



21 Only connected diagrams have to be considered for the computation of A. The disconnected diagrams 
can be resummed into the exponential in e ^ %t( ~ Ho+Hin ^\<^) — e -, ( A /^*|$). 
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in the diagram. We have to sum up all these diagrams to get the right answer. The correction 
to the energy of the state can be schematically written as 



5 A = 5A 

k=2,L=0 



where <5A£ denotes the contribution of diagrams where k fermions are involved in the inter- 
action and L is the number of loops. 



+ 



Figure 2: 2-body, 3-body etc. interactions at tree level and higher loop interactions. 



We will now see how diagrams of different types scale in our limit. Going back to the 
single graviton exchange between a pair of fermions, but now summing over all possible pairs, 
we find that the correction to the energy is 

d 1 

where we used that Np ~ c d +! and m ~ c d +!. In other words 

— — 0(l) 

so the correction is of the same order as the original energy of the system A and remains im- 
portant even in the large c limit. The suppression that we found earlier has been compensated 
by the large number of pairs that we have to sum over. On the other hand, loop diagrams 
between pairs of fermions are suppressed by additional power of \ and are unimportant in 
our scaling limit. 

For a general diagram with L loops, k fermions involved in the interaction (i.e. k ingoing 
and k outgoing ones), Vp fermion-fermion-graviton vertices, and V n vertices of n-graviton 
type, we have the following topological constraints 

V F + Y^nV n = 2I a , 2V F = 2I F + 2k , L = I F + I g - V F - ^ V n + 1 

30 



where I gj Ip denote the number of internal graviton and fermion lines. Now let us count 
powers of c of these diagrams. Each external line contributes a factor of (Nprrif) which in 
our scaling limit is proportional to c. Each n-graviton vertex contributes a power of c 1- ^ and 
each of the Vf vertices a factor of c~5. So in total the powers for a diagram are 



Using the three conditions that we found above this becomes 



So we find that that only L = (tree level) diagrams contribute significantly at large c. The 
conclusion is the following: in our double scaling limit we have contributions from diagrams 
of all orders due to the non-linear nature of general relativity, but only tree level ones. In this 
sense our system becomes (semi) -classical, though nonlinear. 

3.3 From quantum states to fermionic fluid 

Clearly summing up all such diagrams between fermions is not the most efficient way to 
analyze the problem and an approximation method would be desirable. In this section we 
review some basic ideas which are often used in treatments of many-body systems and which 
may clarify the conceptual steps in going from the quantum many-particle description in the 
bulk to the hydrodynamic description and the TOV equations of the next subsection. The 
reader who is not interested in these issues can skip this subsection and go directly to section 



In general, in a many-body quantum problem we first want to determine the ground 
state wavefunction and energy, and then perhaps to study small excitations around it. The 
simplest example to consider is Np non-interacting fermions moving in an external potential. 
In this case we first compute the energy levels for a single fermion in this potential and 
the corresponding wavefunctions ipi{x). The ground state of the many- fermion system is then 
described simply by filling up these orbitals consistently with the Pauli exclusion principle. 
The wavefunction is a Slater determinant of the one-particle wave functions 



5A L k ~ c fc +S(i-f )K 



2 




Iplfa) ^2(^2) ■•• 1pN F {x2) 
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and the total energy is given by the sum of the individual energies 

N F 

E t ot = ej 

When we turn on interactions between the fermions the situation is not as simple, since it 
is no longer true that the many-body state can be constructed from the knowledge of single 
particle energy states. More precisely, single particle energy states, or "orbitals" are not even 
well defined in an interacting system of fermions. By adding or removing a fermion we affect 
the wavefunctions of all other fermions in a complicated way, which in turn influences the 
available energy levels for the fermion under consideration. Hence in principle we have to 
solve the coupled problem of determining the wavefunction for all fermions at once. This also 
means that the ground state cannot necessarily be written as a Slater determinancy. 

A typical problem of this type is finding the ground state of a polyelectronic atom. That 
problem is quite similar to finding the ground state of our fermionic star, where the role of 
the attractive Coulomb force from the nucleus is played by the AdS gravitational potential 
and the electromagnetic interactions of the electrons correspond to the gravitational force 
between the fermions (though the latter is attractive). As we know, solving Schroedinger's 
equation for a polyelectronic atom is not possible analytically and we have to address the 
problem using various approximation methods. We have to do the same for our system. 

In our system the interaction between any two particular fermions is negligible in the large 
c limit due to the - suppression. However since the number of fermions is also very large 
the effect of all of the other fermions on a single fermion can be appreciable. This type of 
many-body problem can usually be treated by mean field theory methods. We replace all 
interactions on a given fermion by an effective external field, which has to be determined in 
a self-consistent way, as we explain below. 

3.3.1 The Hartree-Fock approximation 

In many-body problems the Hartree-Fock approximation is often employed. Its basic assump- 
tion is that the ground state of the many-fermion system can be written as a single Slater 
determinant |\?) of one-fermion orbitals in some appropriately chosen external potential, which 
is determined in a self-consistent manner. One then minimizes the norm of (^>\H\ty) with 
respect to the single-particle orbitals, where H is the full Hamiltonian. This leads to the 
Hartree-Fock equations for the single-particle orbitals ipi and the corresponding energies 6{. 



2 However it can be written as a superposition of Slater determinants. 



32 



In general the physical meaning of the single-particle energies e« in the Hartree-Fock ap- 
proximation is not transparent: for example if we try to add one more fermion to a previously 
unoccupied state €j, then the energy of the new state will not be exactly equal to E tot + €j. 
This is due to the fact that when we add a fermion all other fermions will rearrange themselves 
and in principle one has to solve the Hartree-Fock problem again for Np + 1 fermions, which 
will introduce corrections to the energy, additional to the naive estimate E tot + €j. For the 
same reason exciting a fermion from the state q to the state ej will cost energy equal to €j — 6j 
plus corrections due to rearrangement. 

Nevertheless in the limit of large number of fermions the aforementioned corrections due to 
rearrangement become negligible, at least for certain many-body systems. This is sometimes 
called "Koopmans theorem", which basically states that in the limit of large number of 
particles we can treat any given one of them as a "probe" and define its energy as if it 
had negligible backreaction on the rest of the particles. Hence in this limit the excitation 
energies of a single fermion are simply determined by the Hartree-Fock spectrum e^, by taking 
the difference of the energy of the final minus the initial state. Similarly if we think of 
constructing the ground state by gradually adding fermions and rearranging them to their 
lowest available energy state, then ep would be the energy of the last added fermion. Notice 
that these statements hold only if we talk about changes to a few numbers of fermions relative 
to the total number of them. 

Within the context of the Hartree-Fock approximation we can define a notion of single 
particle density of states, even for a system of interacting fermions 

dN 

9^,^f) = ~^{e,e F ) (3.3) 

which is now a function of the "Fermi energy" €f- This function encodes the density of 
the single-particle wavefunctions of the Hartree-Fock solution with corresponding energies 
6j. We call Fermi energy e F the highest energy of the occupied states and N(e,e F ) is the 
number of occupied single-particle states with energy lower than or equal to e. Clearly we 
have Np = N(eF,ep). In the limit of large number of particles, where Koopmans theorem 
or the probe approximation holds, the function g(e,ep) can be determined operationally by 
doing "spectroscopy" on the state, as follows: we consider the system in its ground state. We 
hit the state by a photon coupled to the fermions and measure the absorption spectrum. This 
gives information about the energy differences between the single fermion states (assuming 
that Koopmans theorem holds) and in this way we can reconstruct the function g(e, €f). 
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The total number of fermions can be written as 

N F (e F ) = f F deg(e,e F ) (3.4) 

where e m in is the lowest of the single-particle statea^j. But now it is important that the total 
energy of the system E to t ^ f F dee g(e,e F ). This is due to the interactions between the 
fermions and the rearrangement issues which cannot be ignored if one tries to build up the 
ground state by adding fermions gradually, since this process involves all of the fermions and 
clearly goes beyond the "probe" limit. However this can still be done for the last added 
fermion, so the following is still true 

dE tot 

— — = e F g(e F ,e F ) (3.5) 
de F 

It is clear that for our problem, that is, for determining the ground state of a large number 
of gravitationally interacting fermions in AdS, the Hartree-Fock method is applicable and 
moreover the assumptions of Koopmans theorem are satisfied. This means that the Hartree- 
Fock orbitals should be really thought as the possible quantum states of a fermion moving in 
the gravitational background produced by rest of the fermions and the density of states (13.31) 
describes the distribution of fermions with various energies, generalizing the non-interacting 
density of state (]2.13|) . 

In practice we will not apply the Hartree-Fock method, since it is still too complicated 
when the number of fermions is large. As we explain below, a further approximation in 
addition to the Hartree-Fock can be made, which leads naturally to the TOV equations. 

3.3.2 Fermion fluid and Thomas-Fermi approximation 

Even after imposing the Hartree-Fock approximation we can make a further simplification as 
follows: instead of solving for the ground state Slater-determinant wavefunction, we can only 
look for the fermion number density n(x) in the ground stated 

n(x) = J d d x 2 ...d d x N \ty(x,x 2 ,...,x NF )\ 2 (3.6) 

which satisfies J d xn(x) = N F . Of course n(x) carries less information than the actual 
wavefunction but it is easier to determine than We express the kinetic and potential energy 



23 Notice that due to the interactions the energy e m i n can be low, for example lower than the rest mass of 
the fermions. 

24 In this section we think of fermions moving in flat d-dimensional space. If we are in curved spacetime 
then the appropriate measure factors have to be included. 
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of the fermions in terms of n(x)o Minimizing the energy with respect to n(x) determines 
the fermion density profile and the energy of the ground state. This is the Thomas-Fermi 
approximation. In many systems the Thomas-Fermi approximation becomes asymptotically 
exact in the limit of large number of particles. As we will explain later, the hydrodynamic 
approximation and the TOV equations in the bulk is the equivalent of the Thomas-Fermi 
approximation (or rather its relativistic generalization) for the quantum interacting fermions 
of the previous sections, also see references PoT - ES] . 

3.4 TOV equations in the bulk 

After these general comments we return to what we would have naively done to deal with the 
gravitational backreaction of a large number of fermions: we revisit our hydrodynamic formal- 
ism of section [231 but now taking into account the gravitational backreaction of the fermionic 
fluid. We will explain how the hydrodynamic approximation is related to the Hartree-Fock 
and Thomas-Fermi approximations in the next subsection 

We have to solve Einstein's equations with cosmological constant coupled to matter 

where the source is a perfect fluid 

T^ u = (p + p)u^u u + pg^ 

and p, p are determined by the flat space equation of state (12.221) for a degenerate fermionic 
fluid. 

As before we look for static, spherically symmetric solutions, so we write the metric in the 
general form 

ds 2 = -A(r) 2 dt 2 + B{rfdr 2 + r 2 d^l 2 d _ v (3.7) 

We parametrize the functions A(r) and B{r) in terms of two new functions M{r) and x{ r ) as 

r 2 " 
- a- - 

„d-2 



25 For example for a system of many (non-relativistic) electrons moving in an external potential V(x) the 
total energy in the Thomas-Fermi approximation would be 

E TF (n) =K J d 3 xn(x) 5 / 3 + J d 3 xV(x)n(x) + y// d 3 xd 3 y "^"^ 
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9, x ( C d M(r) r 2 V , ^ 

B2 (r) = (l-^/- + p) . (3-8) 

where Newton's constant is contained in the coefficient 

16tt G 
d ~ (d-l)V d ^ 

The Fermi gas is static, and hence we take u t = A(r) and all other components of w M equal 
to zero. In terms of M(r) and x{ r ) t ne Einstein equations read 

dM 



dr 



V d - 1 ^(p + P )rB 2 (3.9) 



dx , r C, 
dr 



dp A' 

These constitute the Tolman-Oppenheimer-Volkoff (TOV) equations. 

We choose as boundary condition M(0) = 0, so that M(r) represents the contribution to 
the mass from the energy density inside a ball of radius r, and the total mass is equal to 
M(R). By Birkhoff's theorem the metric outside the star is given by AdS-Schwarzschild. 
Hence, 

\{r) = 0, M(r) = M, for r > R. 
From (I2.30p it follows that the radius R is determined by 

' C ' M ^-(^- (3.10) 



R d-2 (2 \ m J 

Finally the total fermion number can be computed by the integral 



-R 

N F = V d _ 1 I drr d ~ l B(r)n(r) (3.11) 



o 

Let us now see how we can compute the density of single-particle states g(e, ep) from the 
hydrodynamic configuration in the bulk. The energy e of a fermion should be identified with 
the redshifted energy of a fermionic particle in the bulk. At a given point in the bulk we have 
fermions with various values of the local energy ranging from the fermion mass m up to the 
local Fermi energy p(r). So we can write the following formula for the density of states on 
the boundary 

r R 

N(e, e F ) = V^_i / dr B(r) r d ~ l n 



A(r) 
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g(e,e F ) = ^(e,e F ) (3.12) 

where n(e) was defined in ( 12.21(1 . Notice that the statement that E tot ^ J dee g(e, e F ) is related 
to the statement in the bulk that the total mass is not equal to the sum of the redshifted 
masses of the particles. 

In fact, to show that the relation (13. 5p is still valid one needs to use the Einstein equations! 
One has 

dN [ R , . , / dn dB 

e F — = V^ x / drr^ 1 B e F — + —e F n 
de F J \ de F de F 



n \de F de F B 

R dre-(^ + ^ X ') = d -^. (3.13) 
,, \ de F de F J de F 

Here we used the identity (12.231) . the relation (12.281) . the TOV equations (13. 9j) . the definitions 
of A and B in terms of the functions M and \ an d the fact that at the edge of the star 
n(R) = 0. Notice that above the local particle number density depends on the radius as 
n(r) = n (e F /A(r)). 

In the next section we will use the TOV equations to calculate the corrections to the mass, 
the particle number and the density of states defined through the relations ( 12. 12j) . 

3.5 Relation between quantum description and TOV approach 

What is the precise relation between the analysis in terms of quantum states and Feynman 
diagrams in section I3T21 and the hydrodynamic TOV equations (13.91) ? The Einstein equations 
coupled to a perfect fluid are in a certain sense the mean field/Thomas- Fermi approximation 
for a system of fermions interacting gravitationally. This correspondence is of course implicitly 
used in the standard treatment of astrophysical neutron stars, in which the star is treated as 
a classical self-gravitating ball of fluid governed by a fermionic equation of state and not as 
a quantum bound state of individual fermions. The logic in going from the quantum bound 
state to the fluid description was outlined in subsection 13.31 

The equivalence between these two approaches can be made more precise. In [26||27] 
the following problem was analyzed: the authors considered the Schroedinger equation for 
N F fermions mutually interacting with a Newtonian attractive gravitational potential. They 
passed to the Thomas-Fermi description^ replacing the ground state wavefunction by the 



^With a relativistic kinetic term. 
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fermion density n(x). They then showed that in the ground state the fermion density satis- 
fies the (non-relativistic limit of) the TOV equations. In general relativity one should also 
include higher order forces in the Schroedinger equation. Presumably the Thomas-Fermi 
approximation would then reproduce the full relativistic TOV equationscl 

In (28] a related analysis was followed: the starting point is a static spherically symmetric 
metric parametrized by two arbitrary functions A(r),B(r) as in (13. 7p . One solves the Dirac 
equation for fermions moving in this background and places Np of them in the lowest available 
states. Then one computes the stress energy tensor produced by these fermions and uses it as 
a source for Einstein equations for the metric (13. 7p . Then one shows that in the limit of large 
Np solving this problem self-consistently is equivalent to solving Einstein equations coupled 
to a perfect fermionic fluid, that is the TOV equations (13. 9p . 

The reason that we mention all these issues is the following: as explained above, in the 
bulk we have an intuitive understanding of how to start with the microscopic many-fermion 
quantum description and truncate it to the Thomas-Fermi type approximation of the TOV 
equations coupled to a fermionic fluid. What is the equivalent procedure on the boundary? In 
the boundary CFT the description in terms of individual quantum states for the fermions is 
easier to understand, at least in the c — > oo limit before taking the backreaction into account. 
It would be very interesting to perform the equivalent "coarse graining" on the boundary and 
to introduce the boundary analogue of the density of fermions (13 .6p and the TOV equations. 
The main obstacle in this direction is understanding how to deal with the "radial direction" 
in boundary language, see also the discussions in subsection 12.3.51 We hope to revisit these 
issues in future work. 

3.6 Quantum states vs classical geometries in AdS/CFT 

Before we close this section we would like to explain an important conceptual point. In 
AdS/CFT we have an equivalence of two quantum systems, which means that their Hilbert 
spaces are isomorphic. In particular every quantum state on the boundary should be dual 
to a quantum state in the bulk. On the gravity side and at low energies, the Hilbert space 
can be approximated by the Fock space generated by the semiclassical quantization of the 
supergravity fields. At infinite N these fields become free, so we simply have to canonically 
quantize free fields around an AdS background. This semiclassical quantization leads to a 

27 However we should mention that one qualitative feature of the solution, the existence of a critical mass - 
the Chandrasekhar bound, is also visible in the Newtonian gravity approximation. 
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Fock space which is isomorphic to the low-conformal dimension Hilbert space of the boundary 
theory, as we discussed in detail in the previous subsections. 

Our many-particle state is a quantum state: it is a many-particle state constructed by 
acting on the vacuum with many creation operators of the semi-classically quantized free 
fields in AdS. Hence it should not be confused with a classical configuration in the bulk, i.e. 
a classical supergravity solution. The latter should be thought of as a coherent state where 
certain modes of the bulk fields have been coherently excited. In this sense our star differs 
from the "boson stars" (see [21] and references therein) and other solitonic configurations in 
the bulk. Instead, a bosonic analogue of our fermionic star would be, for example, a finite 
temperature bosonic gas like the radiation star considered by Page and Phillips in [32] or the 
"geons" of Wheeler [33] . 

Let us explain this in some more detail: the large N limit is usually thought of as a classical 
limit. Of course no matter how large N is, if the energy of the state under consideration is 
low enough then one would still see the quantization and the discreteness of the spectrum. 
Let us now consider classical supergravity solutions in the bulk. Such solutions have energy 
of order N 2 as can be seen by the overall normalization of the bulk action. However the 
opposite is not true: a state whose energy is of order iV 2 (at large N) does not necessarily 
correspond to a classical supergravity solution. Our quantum states are examples of this kind. 
The difference between the two types of states (coherent vs multi-particle quantum states) is 
that while the total energy is the same (order N 2 ), the distribution of the energy to various 
modes is different: in the coherent states we have few modes excited many times, while in our 
quantum states we have many modes excited few times 3- The latter cannot be described by 
a classical field. More generally, it is important to remember that the fermionic perfect fluid 
is not a "classical field" in the bulk, even thought it is described by a classical density profile. 

Perhaps a useful analogy to keep in mind is the difference between a classical electric field 
and black body radiation. The former has non-vanishing energy density at the classical level. 
By tuning the temperature of the black body radiation appropriately as h — > we can end 
up with a finite classical energy density for the radiation. However there is no sense in which 
the black body radiation can be described by a classical electromagnetic field. The reason is 
that the energy is distributed in too many modes unlike what happens in a coherent state. 



2,N 



We are grateful to S. Minwalla for discussions on these points. 
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4 NUMERICAL SOLUTIONS AND GRAVITATIONAL COLLAPSE 



4.1 Numerical Results 

We will now look for static, spherically symmetric solutions of the TOV equations (13. 9p . These 
equations cannot be solved analytically and we have to resort to numerical integration. We 
present the numerical analysis for AdSs, i.e. d = 4. The qualitative features of our results 
are the same for other values of d. As we explained in section I37L1 we are working in the limit 
G$ £~ 3 — > , Ml — > 00, or in CFT language c — > 00, A — > 00. It is more convenient to present 
the graphs in terms of rescaled dimensionless quantities which stay finite in this limit. For 
the total mass and the total fermion number we define 

c it I cr' b 
and for the fermion mass and local chemical potential in the bulk 

All hatted quantities are kept fixed as c = ^ — >■ 00. 
4.1.1 A typical fermionic star 

We first describe a star of typical values for the total number of fermions Np and fermion 
mass rhf. We plot the radial profile of the local chemical potential in the bulk and the fermion 
number density in figures [3] and H] respectively. The dotted lines in the same graphs show what 
would have been the profiles for the same number of fermions in AdS without gravitational 
backreaction. The edge of the star is the point where (I2.30p is satisfied and the density goes 
to zero. We can see from the graphs that turning on self-gravity moves the fermions towards 
the center of AdS, as expected. 

In figure [5] we plot the function M(r) introduced in (13.81) . which asymptotes to the total 
mass of the solution as r/£ — > 00. We see that turning on gravitational interactions lowers the 
total mass. The difference between the two curves can be understood as the binding energy. 
Remarkably, we see that the binding energy is a very small percentage of the total mass. This 
fact is not a peculiarity of the solution for the the particular choice of Nf,™f but is generally 
true for all of their possible values. In figure |6] we show the #00 component of the metric in 
our solution relative to that of empty AdS where we can see the gravitational well produced 
by the star. 
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Figure 3: Radial profile of the local chemical 
potential with (solid line) and without gravita- 
tional backreaction (dashed line). 
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Figure 5: Mass function with and without 
backreaction. The limiting value as r/i — > oo 
is the total mass M of the system. 
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Figure 4: Fermion number density, with and 
without backreaction. 
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Figure 6: Gravitational potential of the star 
(solid) vs that of empty AdS (dashed). 



4.1.2 Families of stars and a critical mass 

We now want to study a family of stars of increasing fermion number, keeping the fermion 
mass rhf fixed. For this we would like to find the self-gravitating solution of Np fermions in 
AdS, as a function of Np. However when solving the TOV equations, rather than fixing the 
total fermion number Np, it is technically more convenient to fix the local chemical potential 
(1(0) at the center of the star and then integrate the equations outwards. Only after the entire 

(BHD. 
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as a function 



solution has been computed, can the number of fermions be evaluated from 
Following this procedure we plot in figure [7] the mass and fermion numbeic 

of the local chemical potential at the center of the star. 

The most striking feature of the diagram is the existence of a critical value for the mass 

and fermion number which is achieved for a certain value (i c of the local chemical potential 



29 Let us mention that while Nf is of order one in this particular example, it does not mean that we have 
a small number of fermions, since Nf is a rescaled variable. The actual number of fermions is N F — Np c 4 / 5 
which is large for c — > oo. 
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Mo) 

Figure 7: Total mass (blue) and fermion num- 
ber (red) vs chemical potential at the center of 
the star. 
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Figure 9: Binding energy as a function of 
fermion number. 



N F 

Figure 8: Mass vs Fermion number with (blue) 
and without (red) backreaction. 
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Figure 10: Density of states at the Fermi en- 
ergy ip. The value of e F where the density 
diverges corresponds via relation ( 14.1 j) to the 
value of fi(0) where the mass reaches its high- 
est value. 



at the center of the star. This is the Chandrasekhar, or Oppenheimer-Volkoff limit. It will be 
discussed in more detail in the next subsection but for now there are two important points 
that we want to emphasize: first, the right way to read figure [7] is that there are no static, 
spherically symmetric solutions of the TOV equations with a larger number of fermions than 
that achieved by the solution at \x c . If we insist on starting with initial conditions on a 
spacelike slice, corresponding to a larger number of fermions, the solution will inevitably 
become time-dependent and will presumably collapse. Second, as we will explain below the 
solutions with /}(0) > jl c are unstable under radial perturbations and hence in a certain sense 
they should be considered unphysical. We will return to these points in the next subsection. 

In figure [8] we show the total mass as a function of the total number of fermions (keeping 
rhf fixed). In the same graph we can see what would be the total energy M Q of the same 
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number of fermions if they were placed in AdS without taking the self-gravitation into account, 
as computed from (12.101) and (12. lip . The difference between the two curves is very small. In 
figure [9] we see the binding energy relative to Mq. These diagrams reveal a surprising feature 
of our system: the self-gravitating solution breaks down and becomes unstable at a point 
where the binding energy is only a few percents of the total energy. 

Finally we consider the interpretation of our solution in the Hartree-Fock picture on the 
boundary. The Fermi energy on the boundary ep is defined by equation (12.281) so we have 



In figure [10] we plot the density of states at the Fermi energy as a function of the Fermi energy 
£p, which can be computed from (13 .4p . We see that approaching the critical point the density 
of states blows up, which signifies a break-down of the Hartree-Fock approximation. 

4.1.3 Dependence on the fermion mass 

We now consider how the previous results depend on the fermion mass rhf. In figure [UJ we 
plot the total mass as a function of the fermion mass and the chemical potential. We see that 
the critical mass and the value fi c at which criticality is achieved depend on the fermion mass 
rhf. The qualitative behavior is the following: for small fermion mass criticality is achieved 
for higher values of the chemical potential, where the fermions are in the relativistic regime, 
and the radius of the star (which we have not plotted) is large relative to the AdS radius. For 
heavier fermions we have criticality in a regime where the fermions are non-relativistic and 
the size of the star is smaller than the AdS radius. In figures [T2] and [13] we plot the mass and 
density of states respectively, as a function of the Fermi energy and fermion mass. 

4.2 The Chandrasekhar limit 

The most important qualitative feature of our numerical results is the existence of a maximum 
value for the mass of a degenerate fermionic star. This value is achieved for a specific critical 
local chemical potential fi c at the center of the star. It is possible to find static spherically 
symmetric solutions with higher values of density at the center if we take /a(0) > /t c , however 
it is important to notice that increasing /t(0) above fi c does not correspond to increasing the 
total number of particles. As we can see from figure [7] the number of particles has a maximum 
at the same value of /t(0) where the mass reaches its critical value. Solutions with /t(0) > fi c 
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Figure 11: Total mass vs chemical potential Figure 12: Mass vs Fermi energy and fermion 
and fermion mass. number. The underlying pink graph is the re- 

sult without self-gravity. 




Figure 13: Density of states at the Fermi energy vs Fermi energy and fermion mass, plotted 
against the graph without backreaction. 

correspond to a different radial distribution of roughly the same number of particles. Moreover 
as we will argue below solutions with /t(0) > \x c are unstable under radial perturbations. 

What happens if we insist on looking for solutions with a larger number of particles? We 
could do that by starting with initial data on a spacelike hypersurface with a total number of 
fermions higher than the critical one. Then it is clear that these initial data will evolve into a 
time-dependent solution describing a collapsing star. This is the point where the degeneracy 
pressure of the fermions can no longer balance the gravitational attraction, in other words 
the "Chandrasekhar limit" . 

Let us now turn to the stability of solutions with /t(0) > fi c . To check the stability we have 
to linearize the equations of motion around the solution and compute the frequencies of the 
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normal modes oo n , which describe small oscillations around the solution with time dependence 
of the form e~ luJnt . In general the frequencies are real and the perturbations are simply 
oscillating with time, while their amplitude remains bounded. If however there is a linearized 
mode with imaginary frequency, that is with u 2 < 0, then we see that the perturbation will 
grow exponentially with time implying that the initial solution was unstable. 

While we have not performed the linearized analysis, we summarize a general argument 
(see for example [M]) which indicates that the solutions with /t(0) > fi c are unstable. Let 
us say that we have a one-parameter family of static solutions characterized by the central 
density /i(0). If we have a critical point ^ = at /t c , then around that value, to first order, 
we have two different solutions one with fi c and one with fi c + 5fi with the same total mass. 
The change of the profile induced by <5/x(0) can be thought of as a linearized perturbation of 
the solution at fi c . From the fact that the total mass does not change we conclude that for 
this perturbation at fi c we have u> 2 = 0, i.e. we have a zero mode. Generically we can assume 
that oj 2 changes sign as we cross fi c . If we assume that the solutions for /t(0) < fi c are stable 
then we conclude that u 2 < for /}(0) > fi c . Hence these solutions will develop a tachyonic 
mode past the critical point fi c . We expect that at the subsequent critical points more modes 
may become unstable. 

4.3 Endpoint of the collapse 

In a realistic theory the endpoint of the collapse of a degenerate star depends on the details 
of the particle spectrum, interactions, equation of state etc. For example in the real world a 
white dwarf can collapse to a neutron star which in turn may collapse to more (hypothetical) 
exotic states such as quark stars and eventually a black hole. Since we have been mainly 
concerned with generic features of degenerate stars in AdS we cannot answer this question 
in detail. However we notice that, at least for light fermions, the critical mass is of the same 
order as that of a big black hole in AdS. It is reasonable to expect that such an object is the 
most entropic one, for given total mass. This suggests that the final endpoint of the collapse 
will be a big black hole in AdS. 

Nevertheless we should mention that while the endpoint of the gravitational collapse is 
very likely a black hole, the onset of the collapse itself does not have to do with black hole 
physics. The onset of the collapse is marked by the failure of hydrostatic equilibrium for 
a large number of fermions and can be understood in terms of the low-lying (super) gravity 
modes. 
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In principle it would be straightforward to compute numerically the time-dependent solu- 
tion corresponding to the collapse of a fermionic fluid in AdS towards a black hole. In practice 
this problem is significantly harder than what we have done so far, since it involves solving 
partial differential equations instead of ordinary ones (since the fields will now depend on 
time and radius). We postpone the analysis of the time-dependent solutions and the black 
hole formation to future work. 

4.4 Scaling regime 

In this section we will explore some features of the solutions deep inside the unstable regime 
//(0) > fi c . 

4.4.1 Limiting solution 




Figure 14: Profile of the local Fermi momentum kp{r) for increasing values of fi(0) and the 
convergence to a limit curve. 

At very large values of the chemical potential at the center one finds that the solution 
starts to exhibit an interesting behavior: most of the radial profile of the chemical potential 
approaches a limiting shape, while a sharp "spike" develops at the center of the star. This is 
shown in figure [TU where we plot the local Fermi momentum as a function of the radius. The 
Fermi momentum is very large at the center but rapidly drops to the limiting profile. While 
the spike becomes sharper and sharper as we increase /2(0) the total mass contained in it is 
finite (the same is true about the number of fermions in the spike), hence the total mass of 
the solution remains finite. 

One might worry that the spike is an artifact of the numerics, but one can actually find 
an analytic solution which describes its limiting shape. In the limit where fi(0) — > oo we can 
assume that all quantities have a simple scaling behavior as r -> 0. Plugging a power-law 
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ansatz into the TOV equations and looking for a solution in the regime r — > (where the 
cosmological constant can be neglected) one can fix the scaling solution as 

112 77 G r z 



7 

[r) = 
A 2 (r) = r*b 3 + 0(r 



B 2 (r) = - + 0(r) ( 42 ) 



where the constant 63 cannot be fixed by a local analysis^!. One can check that the numerical 
solution of the spike approaches this analytic form in the limit fi(0) — > 00. 

The fact that in the small r limit the TOV equations admit exact scaling solutions in 
which the central density goes to infinity is a well known fact, see for example [35J. 

4.4.2 Oscillations of the mass 

The second interesting feature of the regime /t(0) ^> jl c is the following 31 ! : in figure [7] we 
see that in the limit /t(0) —> 00 the mass M approaches a limiting value M c However 
we notice that the convergence is not monotonic, but rather the function M(/t(0)) undergoes 
small damped oscillations around the critical value M c . In figure [T5] we plot the same function 
at a different scale, in which the oscillations become more visible. This qualitative behavior 
appears in many problems of similar kind. 

In the limit fi(0) — » 00 these oscillations can be described by the following formula 

M{p) w M c + Ae~^ ogm cos(wlog/i(0) + 5) (4.3) 



30 More generally, the scaling solution in d + 1 bulk dimensions is 

p(r) = ^(b 1 + 0(r)), B 2 (r)=b 2 + 0(r) 1 A 2 (r) = r^r (p 3 + 0{r)) 

where 

1 (d-l) 2 (d+l) (d-lf + d-3 

01 = Tl TT^— I o! & 2 = 



AttG (d-l) 3 + d-3' (d-l) 3 -3d+l 

and the constant 63 is arbitrary. 

31 This section was added to our paper after similar calculations were carried out in collaboration with S. 
Bhattacharyya and S. Minwalla in [15]. K.P. would like to thank S. Bhattacharyya and S. Minwalla for very 
useful discussions on these issues during the collaboration in [15]. We would also like to thank V. Hubeny 
and M. Rangamani for very useful comments and suggestions regarding the scaling regime and for bringing 
relevant literature to our attention. 



32 This limiting value depends on the fermion mass rh 



J- 
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log Mo) log Mo) 

Figure 15: Detail of[7]where we see the damped Figure 16: In this plot the oscillations become 
oscillations around M c w 1.1230. This plot more visible, 
corresponds tom/~ 0.95. 

By fitting to the numerical data we find the following values for these constant^ 

M c « 1.230, A w 0.582, 7 « 0.596, w « 0.973, 5 « 3.508 

We will now explain that the formula ( 14. 3 p can be justified by analytic methods and the 
damping constant 7 and frequency u can be determined exactly. Since this aspect is not 
central to the rest of the paper we will be brief and refer the reader to [T5 |[36T438] for more 
details. The main point is to consider the TOV equations as a dynamical system with respect 
to the variable r. In this sense, the critical solution mentioned in the previous subsection acts 
as an attractor fixed point (in the language of dynamical systems) in the space of solutions 
of the TOV equations. The critical solution corresponds to /t(0) — > 00. Solutions with very 
large but finite fi(0) are in a sense small perturbations of the critical solution and hence their 
behavior can be understood by analyzing the TOV equations in a linearized approximation 
around the critical solution. For small r the critical solution is known analytically (14. 2 j) and 
similarly the eigenvalues of small perturbations around it can be analytically determined by 
performing a linearized analysis of the TOV equations around the scaling solution. For our 
system the eigenvalues turn out to be 

A = --±z— (4.4) 

By performing the matching carefully we find that these eigenvalues translate into the follow- 
ing analytic values for the damping constant and frequency of the oscillations of the mass 

3 2V6 
7= 5 ' U = — 



33 The precise values of the constants M c , A, S will in general depend on the value of the fermion mass rhj. 
The reported values are for rht sa 0.95. On the other hand, the constants 7,0; are fixed in an m/-independent 
manner and can be determined analytically as explained in the rest of this section. 



48 



which are in good agreement with the numerical values mentioned above. The small discrep- 
ancies are presumably due to various inaccuracies in our numerics. 

Before we close this section let us mention that while this oscillatory behavior is clearly 
interesting, one should keep in mind that it takes place in the regime /t(0) ^> fi c where the 
solutions are unstable. It would be fascinating if there was a boundary interpretation of the 
scaling behavior and the oscillations of the mass. 

5 BOUNDARY CFT 

In this section we return to the boundary interpretation of the gravitational interactions 
between the fermions. As we discussed in section [2J in the free limit the star is dual to a 
composite multitrace operator (12.2j) . If we first fix the size of the operator and then send c to 
infinity, this correspondence is precise. While the multitrace operator (12. 2p is not protected 
by supersymmetry (i.e. it is not a chiral primary), it is protected due to large c factorization: 
each of its constituents are (descendants of) chiral primaries and are thus protected. Moreover 
at infinite c the conformal dimension of composite operators is simply given by the sum of 
individual dimensions. 

The interactions between the constituents start to become important in the regime where 
we scale the (bare) dimension of the operator as Ao = ec and keep e finite. Now we want 
to understand these interactions from the boundary point of view. We want to work in a 
regime where the boundary theory has a semiclassical gravity dual, which implies that the 
gauge coupling must be large (for example large A 3> 1 in the Af = 4 SYM). This means that 
we cannot analyze the boundary theory perturbatively in the gauge coupling. However even 
at strong 't Hooft coupling the theory has an expansion in -. As we will see the =■ expansion 
together with certain basic assumptions about the boundary CFT can be used to capture 
some of the qualitative features of the system that we found in section HJ 

In the free theory the conformal dimension of the operator (12. 2ft is dual to the energy of 
the fermionic gas in the bulk. The correction to the energy of the gas due to gravitational 
interactions, that is the binding energy, is related to the anomalous dimensions of operators 
of the form (|2.2p . once - corrections are taken into account. When e — > these corrections 
become negligible but they are important when e is of order 1. In a certain sense these 
corrections are controllable as we turn on e slowly. 

The precise form of the - corrections to the composite operators will depend on the details 
of the conformal field theory. So far we have mainly considered the neutron star in AdS, in 
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a model where the only degrees of freedom consist of a massive fermion coupled to general 
relativity. Of course, one of our main interests was to determine the dual description of the 
neutron star and in particular to understand some microscopies of the neutron star instability 
and its presumed collapse into a black hole. In order to achieve this, we need to embed 
our discussion in a proper AdS/CFT duality, like e.g. the duality between Af = 4 SYM 
theory and type IIB string theory on AdS 5 xS 5 . However, one immediately see that this is 
potentially difficult. The mass of the fermion £A scales like N 2 / 5 £ in the limit we consider, 
which is larger than the typical mass of an excited string state 1/4 ~ (A) 1 / 4 ^ in the 't Hooft 
limit. Therefore, one expects that both massive string degrees of freedom, as well as a large 
number of Kaluza-Klein modes of the massless string degrees of freedom, should be included 
to properly discuss the physics of the neutron star. 

Before discussing these issues in more detail, it is instructive to imagine a situation where 
these additional degrees of freedom have somehow been decoupled. Though perhaps unreal- 
istic it still interesting to see how much physics such a toy model could capture. In the rest 
of this section we will focus on a "toy model" CFT, where the only relevant fields are the 
graviton and the fermion, up to a cutoff of the order of the Planck scale in the bulk. 

5.1 Graviton exchange and first correction to anomalous dimen- 
sions 

In the toy model setup the bulk physics is completely described by a massive fermion coupled 
to general relativity. The corresponding field theory statement is that the only low-lying 
"single trace operators" in the field theory are the operator ^ and the stress tensor T^, 
together with their conformal descendants. Of course we also have to include their multi- 
trace composites, as required by crossing symmetry. However these do not correspond to new 
fields in the bulk, but rather to multi-particle excitations of the basic fields. For simplicity 
we will restrict attention to the case d = 4 only. 

In this simple context, it is possible to estimate the point where the instability will ap- 
pear, using known properties of the operator expansion of conformal field theories in four 
dimensions. The idea is to consider the anomalous dimension A of the composite operator 
$ as a function of particle number Np, and to determine when dA/dNp changes sign, which 
is the direct analogue of the computation we did in the bulk. Now in order to compute the 
anomalous dimension of the operator we are going to use the fact that it is built out of 
operators \1/ and their derivatives, and the relevant anomalous dimensions can be extracted 
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from correlations functions of \I> with itself. To compute these correlation functions, we will 
use the fact that in the OPE of \I/ with itself only the stress-tensor can appear to leading 
order in -, apart from multi-particle composites made out of This simplification may not 
be accurate in a complete unitary CFT but by assumption this will happen in our toy model. 
Altogether this computation will therefore be a direct field theory analogue of graviton ex- 
change, where instead of graviton exchange we have the exchange of a stress tensor in the 
intermediate channel of an OPE. 

For future purpose we discuss the relevant computations in a slightly more general setting. 
Thus, we consider two self-adjoint conformal primaries Oi,0 2 . Using the OPE we can define 
a composite operator : 0\0 2 '■ whose conformal dimension is Ai + A2 to leading order in 1/c 

1 (x)0 2 (y) = ...+ :O l 2 :(y) + ... (5.1) 

This form of the OPE, and the presence of an operator with the quantum numbers and 
dimension of : 0\0 2 : can be extracted by performing a conformal block decomposition of the 
4-point function, if we assume that at infinite c it factorizes to a product of 2-point functions. 
At the first non-trivial order in 1/c the dimension of the operator : 0\0 2 '■ will be 

A 12 = A! + A 2 + ^ (5.2) 

c 

where 5i 2 is the "anomalous dimension". To extract 5\ 2 , we consider the 4-point function 

(0 1 (x 1 )0 1 (x 2 )0 2 (x 3 )0 2 (x A )) = * + 0(1) (5.3) 

F12I F34T c 

Here we assumed that the two-point function of 0% and 2 vanished to leading order in 1/c. 
Now let us see what happens at order -. As explained above, we will only consider the effect 
of a stress-tensor in some intermediate channel, e.g. the (12) — > (34) channel. The appearance 
of T in the OPE of O1O1 is controlled by the Ward identities and we have schematically [39] 



Oi{x)OM ~ ... + Ax + ... (5.4) 



while 



and also 



2 (x)0 2 (y) ~ ... + A 2 T + (5 5) 

i x y 1 



T(x)T(y) ~ j-H-p (5.6) 



51 



Oi(xi) 



2 (x 3 ) 



T, 



1 (x 2 ) 



2 (x 4 ) 



Figure 17: Conformal partial wave correspond- 
ing to the exchange of T^ v in the double OPE. 
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Figure 18: Graviton exchange Witten diagram. 



From this we can compute the contribution of T and all of its descendants in the double 
OPE (12) — > (34) and we can write the result in terms of a "conformal partial wave." The 
final result is given in eqn (4.3) in [39|. Taking this result, expanding it in the (13) — » (24) 
channel, and isolating the logarithmic singularity in that channel leads to 

AxA 2 



'12 



4c 



For e.g. N = 4 SYM we have in the conventions of [39] c 
is 

#12 ' 



(5.7) 



so the anomalous dimension 



AiA 2 



N 2 



(5.8) 



This has exactly the same dependence on Ai, A 2 and N as one would expect from a bulk 
graviton exchange between two particles of mass £A\ and £A 2 separated a distance r ~ £ i.e. 



G 



ni\m 2 



This strongly suggests that one should also be able to determine that to leading order in 1/c 
the conformal dimension of $ behaves as 



A($) 



a A 2 

A - a W 2 



(5.9) 



where a is some number of order unity. The result (15.91) would support the bulk result 
that the neutron star becomes unstable at A ~ iV 2 . It is however not straightforward to 
establish ( 15. 9ft . One first needs to generalize the above computations to correlations functions 
involving ^'s and their derivatives, and the latter are no longer conformal primaries. In 
addition, in these four-point functions not only T but also bilinears in \1/ and its derivatives 
will appear in the intermediate channel, so that additional conformal partial waves need to be 
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included in the computation. These additional bilinear contributions to the double OPE are 
necessary in order to have a crossing-symmetric 4-point function at order the conformal 
partial wave corresponding to exchange is not crossing symmetric by itself. It has to be 
"dressed up" with contributions from the exchange of fermion bilinears in order to combine 
into a crossing-symmetric contribution which basically becomes the graviton exchange Witten 
diagram. Finally it would be important to clarify possible complications from the mixing of 
the composite operators with other operators of the same quantum numbers, once interactions 
are taken into account. 

Apart from these complications it seems likely that the single graviton exchange in the 
bulk can be reproduced by the T exchange in the double OPE. As we mentioned, the latter is 
fixed by conformal invariance and the Ward identities, so it can be reliably computed even at 
strong 't Hooft coupling if we know the conformal dimensions A 1; A 2 and the central charge 
c. Thus the gravitational force between two particles in the bulk can be understood on the 
boundary via the T exchange, which is a universal property of CFTs. To compute the full 
backreaction in our scaling limit we need to sum over all pairs of particles, but also to include 
the effects of 3-body, 4-body etc. interactions in the bulk as we discussed in section 13.21 It is 
not clear to what extent these higher order forces can be explained on the boundary by some 
universal argument, as was the case for 2-body forces. 

On the other hand, if we have established that the 2-particle force in the bulk is reproduced 
by T exchange, then to compute the total backreaction from 2-particle forces we simply have 
to sum over all pairs (i.e. all pairs of single trace constituents of the multi-trace operator <fr). 
This summation is an involved combinatoric problem but it has nothing to do with strong 
coupling effects in the gauge theory, so in principle it should be tractable. This part of the 
backreaction is the "Newtonian" limit of general relativity, in the sense that we ignore self- 
interactions of the gravitons. While this is not the full answer, we would like to emphasize 
that the existence of a maximum mass and the Chandrasekhar limit can also be estimated by 
treating gravity as Newtonian. 

5.2 Boundary description of the instability and the collapse 

The degenerate star represents a state in the conformal field theory with a large number of 
fermionic "glueballs", placed in a configuration with lowest possible energy. It is a high- 
density state at zero temperature. If we keep increasing the density by adding more glueballs, 
then at some point the state becomes unstable and presumably undergoes a phase transition 
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towards a deconfined "quark-gluon plasma" thermal state, the dual of a black hole in the 
bulk. 

It would be interesting to further explore the physics of this instability and the meaning 
of the tachyonic mode in the CFT. We expect that when the number of single trace operators 
multiplied together becomes very large the resultant operator mixes very strongly with other 
operators in the CFT and eventually with the operators dual to the black hole microstates, 
which have a very large entropy. Presumably the time scale set by the condensation of the 
tachyonic mode contains information about the thermalization time of the dual field theory. 

6 VALIDITY AND EMBEDDING IN STRING THEORY 

Our discussion so far has been general, without reference to a specific AdS/CFT setup. In this 
section we would like to discuss a couple of issues related to the validity of our approximations 
and the possibility of embedding our story in a precise holographic duality. 

There are various possible sources of corrections to our analysis. We will address them in 
order of complexity as follows: first we will assume that the theory in the bulk consists only 
of a massive fermion and the graviton and check the validity of our approximations. Then we 
will include the effect of other massless fields, Kaluza-Klein modes on the sphere as well as 
stringy degrees of freedom. 

6.1 Toy model with a single fermionic field and the graviton 

A possible complication that we have not dealt with so far is that without a fermion con- 
servation law, the number of fermions does not have to be constant in time as there may be 
processes under which fermions convert to gravitons and vice versa. Without a conservation 
law it is natural to expect that if we wait long enough then part of the fermions will be con- 
verted to gravitons and the final state will be a finite temperature thermal gas of gravitons 
and fermions in equilibrium. 

Nevertheless we will now argue that the timescale for this thermalization can be made 
parametrically large, in such a way that the degenerate fermionic star is a good description of 
the system for quite a long timescale. This is based on the ^ suppression of the interactions 
which are responsible for changing the number of fermions and gravitons. Let us now explain 
this point in more detail. We want to estimate the decay rate of the fermion gas to a mixed gas 
of fermions and gravitons. Notice that provided that the typical wavelength of the fermions 
is small (i.e. /t(0) is large), which is true in our limit, this question can be studied locally by 
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looking at a small box in AdS containing the fermionic gas and asking how quickly it converts 
to a mixture of fermions and gravitons. In other words it is a question about the validity of 
our equation of state, which can be analyzed even in flat space. So we consider a box of the 
gas at chemical potential \i and we try to estimate its "lifetime" r. The relevant scales in the 
problem are the fermion mass rrif, the chemical potential fi and the planck mass m p which 
controls the interactions. For simplicity we will only check the validity in the regime where 
fi » rrif, in which the particles are extremely relativistic, and thus interactions are more 
energetic. If the approximation turns out to be reliable in this limit then we expect the same 
to be true for less energetic particles i.e. when fi w rrif. 

In this limit we have only two dimensionful parameters \x and m p . We write the lifetime 
of the fluid as 

\m p J /i 

for some function f(x) of the dimensionless ratio involved in the problem. Our scaling limit 
corresponds to ji ~ c~ while in general m p ~ c^r. In other words x — Y 0. Let us assume 
that / has a smooth limit in that regime and expand 

f(x) = x a + subleading 

How can we determine the exponent a? This limit can also be understood as one where we 
keep m p fixed and send \x to zero. This is a limit of an infinitely dilute gas, where we can 
use Feynman diagrams to estimate the lifetime. As a typical example consider the leading 
Feynman diagram of two fermions going into two gravitons. From such diagrams it is easy to 
determine the dependence of the lifetime on Newton's constant and we find that to leading 
order r goes like ^, which implies a = —d + 1, or the lifetime 

r ~ — — h subleading, for > 

jj, d m p 

or 

tI ~ c*+i + subleading, for c — > oo 

Hence the star becomes long-lived (in AdS radius units) at large c. Similar scaling holds for 
other possible diagrams^, if we assume that the correlators scale with c in the way expected 
from a theory with a standard large c expansion^. 



34 Notice that in the regime where ^> m/ the particles are effectively massless, so from energy momentum 
conservation there are no diagrams with 2 incoming and 1 outgoing particles. 

35 This means that the connected correlator of n "single trace" operators scales like c~ , in a normalization 
where the 2-point functions of single trace operators are of order one. 
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6.2 Including other modes 

The effects of these new fields can be roughly divided in two categories: first they may lead to 
a modification of the equation of state for the fermions (or even to a more drastic breakdown 
of the approximation of a fermionic fluid in the bulk) and second they may introduce new long 
range forces (for example electromagnetic ones if the fermions carry R-charge). The effect 
of new long range forces is easier to deal with: one has to solve the TOV equations coupled 
to additional massless modes. For example we will discuss in more detail the inclusion of an 
electric field in the next section. As we will see the qualitative behavior of our system remains 
the same. In general the number of massless modes which can be sourced by the fermions 
will be small and we expect that generically the gravitational coupling will be the dominant 
ona 36 l. So we do not expect that the inclusion of new long range forces will modify our results 
qualitatively. 

On the other hand the corrections to the equation of state may be more drastic. To be 
concrete, we will first discuss whether these corrections are under control in the usual IIB 
AdSsX S 5 background dual to M = 4 SYM. In this case we run into a difficulty The mass of 
the fermion £Aq and the chemical potential /i scale like N 2 ^ 5 £ in the limit we consider, which 
is parametrically larger than the typical mass of an excited string state l/£ s ~ (ggNy^i at 
weak string coupling. Therefore, we expect that massive string states will be produced by the 
interactions between the fermions. Notice that our scaling arguments of the previous section 
do not apply since we can also have diagrams of two fermions producing a single massive 
string state, for which the 4 suppression is not sufficiently large. Moreover the Hagedorn 
growth in the density of string states (and thus in the density of possible final states) suggests 
that the longevity of the star cannot be guaranteed. Instead the star may quickly convert 
into a very complicated mess of stringy modest 7 .. 

An alternative AdS/CFT realization of our star would be in M-theory backgrounds, i.e in 
AdS4X S 7 or AdSyX S . The main advantage of these theories is that all new physics beyond 
supergravity is related to the single UV scale mn, the eleven-dimensional Planck mass. Since 
we are working in a limit where all energy scales are parametrically smaller than m n we can 
safely ignore all M-theoretic objects beyond those of eleven-dimensional supergravity. To see 
this more precisely, let us from the start consider the full eleven- dimensional picture. In this 



36 One way to explain this expectation is that other long range forces will only couple to the "charges" of 
the basic fermion \Jj, while gravity couples to the mass of "J as well as the kinetic energy of the fermions i.e. 
the "derivatives" acting on 

37 We would like to thank S. Minwalla for very helpful comments about the validity of our approximations. 
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case the massive fermion corresponds to a certain Kaluza-Klein mode of the supergraviton 
multiplet on the sphere. Now let us try some specific expressions. For simplicity we ignore 
factors of order 1 in the rest of this section. In the case of M-theory on AdS4 x S 7 the AdS 
radius £ is related to the eleven-dimensional Planck mass mn by 

£ = - — 6.1 

mn 

where N is the number of M2 branes, or the units of 4-form flux. In other words the central 
charge is c = iV 3//2 . It is not hard to see that in order to have nontrivial gravitational 
backreaction then we need an energy density (from the 11 dimensional point of view) which 
scales like 

So the lid chemical potential of our fluid is of the order 
Plugging into the formula for the lifetime for d = 10 we find 

so it increases with iV (though very slowly). This shows that if we have no new degrees of 
freedom at scales parametrically smaller than mil then in the large N limit our fluid becomes 
stable. 

Let us now explain why we do not expect any such degrees of freedom in M-theory on 
AdS4X S 7 . Consider the sequence of compactifications of M theory on AdS4X S 7 labeled by 
the integer N. In the large N limit the two scales that we are familiar with are the AdS scale 
I and the lid planck mass mn which are related by (I6.ip . Let us now assume that there 
exists another mass scale, which lies somewhere between l/£ and mn where new massive 
degrees of freedom appear, in analogy with the string scale in IIB. Let us call this scale m new . 
We now consider how this scale behaves, relative to the AdS scale, in the large N limit. In 
other words, what is the limit of the product £ m new at large N7 According to the AdS / CFT 
duality for this system we have 

£ TTlnew ^ 

otherwise the statement that in the large N limit we can approximate M-theory on AdS4X 
S 7 by lid supergravity would not make sense. 
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So assuming this condition, we conclude that m new will necessarily have to be a new "UV 
scale" (i.e. a small length scale relative to the size of AdS). So we learn that in the large 
N limit we have two UV scales m new and m\\. Notice that as we make N larger and larger 
both of these scales move even more towards the UV, and become more and more tiny as 
length scales compared to the AdS scale. If in the large N limit we had two such UV scales 
in M-theory on AdS it is natural to assume that the same would be true in flat space, since 
it is hard to imagine how the finite radius of AdS would introduce a new UV scale. 

On the other hand we know that lid flat space M-theory has only one scale mu, so m new 
has to be proportional to mu in the large N limit and not parametrically separated from it. 
This is because there cannot be a free parameter to tune the ratio between the two as there 
are no moduli in lid flat space M-theory. So essentially there is only one UV scale m n in 
the large N limit. Then we can apply the arguments of the previous sections as before. 

This argument suggests that our neutron star can probably be reliably embedded in M- 
theory on AdS4xS 7 . 



In this section we analyze the degenerate star in the case that the fermions are charged under 
a £7(1) gauge field. One motivation to do so is the following: in most known examples of 
AdS / CFT dualities the light fields in the bulk correspond to chiral primaries on the boundary 
i.e. they have nontrivial R-charge. Hence the degenerate star will source the corresponding 
electromagnetic field, which in turn will induce an additional electromagnetic force on the 
fermions. We would like to understand to what extent the inclusion of these forces will modify 
our results. Another motivation is the recent analysis of bulk systems at finite charge density, 



usually in the presence of a black hole, in the context of "holographic super conduct ivityT 8 !. 

How can we determine the equilibrium configuration of a large number of charged fermions 
which are interacting gravitationally and electromagnetically? We would like to work directly 
in the hydrodynamic description and derive the analogue of the TOV equations for charged 
matter. First let us understand what is the equation of state that we have to use. Does the 
addition of charge modify the equation of state (I2.22p . (I2.23p . (I2.24p ? An equation of state 

38 While this draft was being prepared the work [21] appeared, which has some overlap with the results 
of this section. That work focuses on bulk configurations of a charged fermionic fluid in the Poincare patch 
with planar sections (in contrast our stars are in global AdS and have spherical symmetry). Presumably the 
solutions of [3T] can be recovered by taking the planar limit of the solutions constructed here, while scaling the 
fermion mass and charge appropriately at the same time. We would like to thank S. Hartnoll for discussions 
about their work. 



7 CHARGED STARS 
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is defined as the relation between various intensive quantities in the thermodynamic (infinite 
volume) limit. However when we have long-range forces sourced by the fluid it is impossible to 
take the infinite volume limit (the system suffers from the analogue of the Jeans instability). 

To see what is the right way to include the effect of the charge let us remember how we 
incorporated the gravitational interactions of the fermions: we used the standard equation 
of state, but instead of simply solving the hydrodynamic equations V^T^ = on a fixed 
background, we solved the Einstein equations using the fluid as a source, which led to the 
TOV equations. Similarly in the case of a charged fluid what we have to do is to solve the 
Maxwell-Einstein equations, where the fluid sources both the Einstein equations by its stress 
energy tensor and the Gauss law for the U(l) gauge field by its charge density. 

7.1 Charged TOV equations 

Here we derive the TOV equations for charged degenerate stars. We have the Einstein equa- 
tions 

R»u - -jRg^ + = 

where the total energy-momentum tensor is now that of a perfect Fermi fluid plus that of the 
electromagnetic field 



= (p + vWu v + pg^ + (f w F£ - ^g^F^F^j 



(7.1) 



Gauss's law reads 

d^y/gFn = tey/gJ" (7.2) 

Turning on only the radial electric field E(r), we have F i = —F w = ^/g tt g rr E(r) as the only 
non-vanishing components of the electromagnetic field strength. Inserting this T^ u into the 
Einstein equations merely has the effect of shifting the energy density and pressure by the 
following 

E(r) 2 , E(r) 2 
P + — and p -> p — 

87T 87T 

This implies that the form of the first of the equations in f)3.9p gets modified to 

whereas the other two equations in (I3.9P are only implicitly affected through M(r), which 
enters the functions A(r) and B(r) in (13. 8p . Now, conservation of T^ u yields 

t = - {p + p) W) + ^ [{d ~ 2) E{r) + rE ' {r)] (7 ' 4) 
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The electric field enters the system of equations as a new variable, hence we need another 
equation in order to uniquely specify the solution. This comes in the form of Gauss' law, 
written as 

E'{r) = - ( ^" 2)j5;(r) + p ch B(r) (7.5) 
r 

where the charge density p c h{ r ) — <?/ n (?") for fermion charge qf and the number density n(r) 
can be obtained from the thermodynamic identity (12.231) . All the primes denote derivatives 
with respect to r. Finally substituting eq.f l7.5p into eq.f l7.4p gives us the charged analog of 
the last equation in (13. 9p . 

d P ( , \ A '^ , Vd ~ 2 \n( \ (7 r\ 

Tr = - {P + P) M^ + ^- Pch E{T) B{r) (7 - 6) 

The boundary conditions remain the same as in section I3~4l for the uncharged case. Also eqs. 
(13.1 Op and (13.111) for the radius R and total fermion number Np remain unchanged in form, 
except that the functions M(r) and B(r) therein are now replaced by their charged counter- 
parts. We now proceed to solve this system of coupled differential equations numerically. 

7.2 Numerics of charged solutions 

Extending our results of the last section, we now look for numerical solutions of spherically 
symmetric and static, but charged TOV equations for degenerate stars in AdS 5 . As discussed 
in section H~T1 we are working with parameters scaled as G^/i 3 — > and Ml — > oo. Therefore 
other physical quantities get the following rescaling 

where the hatted quantities are kept fixed as c — > oo. Notice that in this normalization the 

1 /2 

charge includes a factor of the U(l) coupling constant which goes like G x 5 '\ This means that 
in the normalization where the charge is integral it scales like c 1//5 which is consistent with 
the scaling c 1//5 of the fermion mass rrif for the case where the fermions are (descendants of) 
chiral primaries. 

7.2.1 A charged degenerate star 

In this section we present solutions of typical degenerate stars comprised of charged fermions. 
In order to study the effects due to varying charge, we keep both the total fermion number 
Np as well as the fermion mass rhf fixed. Then naturally, upon varying the charge of the 
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star, we expect to see a variation in the total mass as well as the chemical potential profile. 
Figure [19] shows the radial mass distribution for four stars with different charges. The mass 
profiles approach the total mass as rjl — > oo. Notice that turning on fermion charge gradually 
increases the mass even though the total number of fermions is held fixed. This is what we 
expect from eq.( l7.3[) . where the square of the electric field enters the mass integral. What 
happens here is that electrostatic repulsion acts against gravitational attraction and that 
reduces the binding energy compared to the uncharged case. 



M(r) 




k F (r)/rhf 




r/£ 

Figure 19: Radial mass distributions of 
charged degenerate stars. Blue profile denotes 
the uncharged case; red, black, orange denote 
profiles with succesively increasing charges. 
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Figure 20: Radial functions of Fermi momen- 
tum in ascending order of charge, with un- 
charged case shown in blue. 
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Figure 21: Electric fields inside and outside 
the stellar mass, scaled by r 3 . Flat line (blue) 
denotes uncharged case, orange curve repre- 
sents maximally charged case. 




Figure 22: Gravitational potential of the 
charged degenerate stars. 



In figure [20] we plot radial profiles of the local Fermi momentum in the bulk for the same 
set of charges as in figure [191 The radius of the star can be read-off from the x-intercept, 
where the pressure, energy and number densities all vanish. Compared to the uncharged star 
(denoted by the blue curve), we see that increasing fermion charge has the effect of increasing 
the overall size of the star. Also the fermionic momentum at the core of charged stars is less 
than that of the uncharged one. Both these effects can be understood as the electrostatic 
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repulsion acting against the gravitational attraction. Remarkably, as we see from these plots, 
this cancellation effect is only small, given that we have determined solutions for a wide range 
of fermionic charges. 

Figure I2T1 gives the solution of the electric field of the same stars, scaled by r 3 . The blue 
profile is zero everywhere, as it should. For the other profiles, we note that they exactly 
flatten outside the radius of the respective stars. For a spherically symmetric metric in five 
dimensions, this means that the electric field outside the star behaves like one resulting from a 
point charge that falls off as 1/r 3 . Finally in figure I2"2l we see the goo component of the metric. 
As expected, the minima of the charged solutions lie slightly above that of the uncharged. 

What we notice from the results presented in this section is that deviations from the 
uncharged star upon tuning on fermion charge are relatively only a small fraction and com- 
parable to the difference between the self-gravitating and unbackreacted star in the previous 
section. 

7.2.2 Family of charged TOV solutions and critical mass 

We now study the characteristics of a family of charged solutions. This will help us identify 
stable solutions for static, spherically symmetric charged degenerate stars. Furthermore this 
analysis also reveals some interesting features of the solution into the unstable phase (after 
the critical point). The results are presented in figure l23l Here we have four profiles in order 
of increasing fermion charges; violet being the uncharged case and green being the maximum 
charge. The fermion mass rhj has been held fixed throughout. Each point on these profiles 
represents a star with a given central chemical potential /t(0). Knowing this, one can easily 
compute the energy density at the core of a star p(0), which is plotted on the a;-axis with a 
logarithmic scaling. 

On the vertical axis of figure [23] we have the mass. Just as in figure [7] in section 14.1.21 
once again we clearly observe the existence of a critical mass, seen at the overall maximum 
of each profile. This is the Chandrasekhar or Oppenheimer-Volkoff limit for charged stars. 
After this limit, the solutions are unstable against radial perturbations for the same reasons 
as discussed in section 14.21 Compared to the uncharged case, we see that the critical mass 
rapidly increases upon turning on a charge. The figure also shows that the critical chemical 
potential fi c gradually increases and so does the amplitude of oscillations in the unstable 
region like those discussed in section 14.4.21 In conclusion, we observe that many qualitative 
features regarding the onset of gravitational collapse discussed in section l4~2l remain valid also 
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for charged degenerate stars. 

These observations also suggest that the holographically dual CFT picture for charged 
stars is a relatively straightforward extension of the uncharged case. In addition to the 
exchange between single-trace constituents we also have to include the exchange of the 
conserved current dual to the U(l) gauge field in the bulk. 

Finally let us also note that the charge of the star cannot be increased indefinitely, keeping 
other parameters fixed. We have checked numerically, that in that event the solution itself 
breaks down. This can be understood in the following way. From the discussion in section 
17.2.11 above we have seen that the electrostatic force of the star acts against its self-gravity. 
Once this repulsion becomes large enough to overtake self-gravity, it inhibits the formation 
of a star. On the other hand, there might be a double scaling limit which would allow us 
to construct arbitrarily large charged stars and to take the "Poincare limit" of global AdS, 
which would most likely connect to the solutions considered in [21]. It would be interesting 
to explore this further. 
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Figure 23: Total mass vs central energy density (logarithmically scaled) for a family of stars 
with increasing values of fermion charge, from uncharged case (violet) to maximally charged 
(green) . 

7.2.3 Attractor fixed points 

Another feature of the limiting solutions deep into the unstable regime is the existence of 
attractor fixed points in the solution space of TOV equations is shown in figure [2H where 
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the mass has been plotted versus the radius of the star for a family of stars each with a 
specific central chemical potential fi(0). If we interpret this in a dynamical systems language, 
the radial direction plays the role of the time parameter. In the scaling limit £i(0) — > oo, 
the solutions remarkably spiral towards the attractor fixed points of the TOV equations. 
In figure [2H the different profiles denote solutions with different fermion charges, with the 
violet curve being the uncharged case and the green profile being maximally charged. These 
spirals are precisely the mass (and similarly radius) oscillations that we discuss in section 
14.4.21 below. Here we see note that these spiraling attractors/mass oscillations are a generic 
feature of the TOV system of equations, both with and without U(l) charges. Figure 1251 
shows a blow-up near one of the fixed points. In these plots the fermion mass rhf was held 
fixed. Upon reversing the set-up to vary rhf and keep the fermion charge fixed, we have also 
observed similar spiraling attractors. In this sense, the fermion mass and charge fully specify 
the parameter space of all fixed points. 




R/£ R/£ 

Figure 24: Mass vs radius of a family of stars Figure 25: Zooming-in the mass vs stellar ra- 

for increasing fermion charges from blue curve dius profile near the fixed point. The blue 

(uncharged) to green (maximum charge). curve denotes a family of uncharged stars. 



8 THERMAL GRAVITON STAR 

In this section we will study the thermal version of a star in an AdSd+ixS fc background. We 
will consider a thermal gas ofalH supergravity modes at sufficiently high temperature so that 
they backreact to the geometry. Our solution will be the analogue of the radiation star [32] , 
including certain modifications due to the internal sphere S fc . 

The main motivation for considering the thermal star is the following: if we start with a 



39 In this section when we refer to a "graviton star" we refer to a system where all massless supergravity 
fields have been turned on thermally. 
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generic, sufficiently complicated, initial state in AdS we expect that at late times 40 ! the system 
will be approximately described by a thermal density matrix of an ensemble with the same 
values of conserved charges (such as energy or R-charge) as the initial state. This final thermal 
state may be a black hole in AdS, or a thermal gas of all supergravity fields depending on the 
initial state. We want to study the thermal gas endpoint in the case where the temperature 
is high enough for the gas to backreact. 

Moreover, the thermal star may describe a "superheated" phase of certain gauge theories 
as we now explain. Large N gauge theories with classical gravity duals undergo deconfinement 
phase transitions [40] , which are the analogue of the Hawking-Page phase transition between 
a gas of gravitons at low temperature and a black hole at high temperature. This phase 
transition takes place at a temperature Tup which is of order one (i.e does not scale with N). 
The low temperature phase has energy of order N° while the high temperature has energy of 
order iV 2 . The jump of energy between the two phases indicates that the phase transition is 
of first order. What would happen if we considered the same system in the microcanonical 
ensemble? Notice that the thermal gas phase, while thermodynamically subdominant, is 
locally stable even at T > Tup- This means that if we work in the microcanonical ensemble 
we can push the thermal graviton gas into a superheated phase by pumping in energy so that 
the effective temperature T of the gas will go above the Hawking-Page temperature ThplL- 

Let us consider what happens to this superheated phase as we increase the temperature T. 
From general intuition we expect that the temperature cannot be increased indefinitely since 
the graviton gas will undergo gravitational collapse, if sufficiently heavy. The energy of a 
thermal gas in AdS^+i grows like T d+l so it will backreact to the geometry when GT d+1 £ 2 ~ f 
where G is Newton's constant. In terms of the central charge we find that gravitational 
backreaction of the thermal gas will start to take place when T ~ c 1 /( ci + 1 )^- 1 i n the large c 
limit] 42 !. So in theories where there is no other scale up to the Planck mass (for example in 
theories without a "string scale"! 43 !) we expect that the superheated thermal gas phase will 
indeed persist up to temperatures of the order where gravitational backreaction will become 
important. We want to compute the effect of the backreaction and understand thermodynamic 



40 Dcpending on the choice of initial state the thermalization may be slow in units of the AdS scale, in the 
large N limit, as we argued in the previous section. However if we wait long enough the system will eventually 
thcrmalize. 

Of course this would not make sense in the canonical ensemble. 

42 Notice that this is parametrically smaller than the Planck scale in the bulk which grows like mp ~ c 3 ^ 1 ^ -1 . 

43 As we discussed in the previous section, this phase is not relevant for the J\f = 4 SYM, since the re- 
quired temperature is much hotter than the Hagedorn temperature but it might be realizable in M-theory 
backgrounds. 
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properties of the superheated phase, for example its equation of state E(T). 

To compute the backreaction of a thermal gas we have to solve the Einstein equations 
coupled to the stress tensor of a fluid at temperature T. For massless fields in the bulk the 
relevant equation of state is that of radiation. If the massless fields have spin we simply have 
to multiply the energy and pressure by the appropriate number of polarizations. This means 
that a thermal graviton star in AdS^+i is the same as the radiation star of [32J, apart from 
a difference in a numerical coefficient in the equation of state, relating the energy density to 
the temperature, due to the number of graviton polarizations. Then we do indeed see that 
there is a maximum temperature beyond which the superheated phase ceases to exist. 

In this section we would like to repeat the computation of [32J in backgrounds of the form 
AdSd + ixS fc i.e with an internal sphere, since they are more relevant for AdS/CFT. For this 
purpose we start directly with the d + k + 1 dimensional supergravity equations, sourced by a 
d + k + 1 dimensional radiation fluid. After describing the structure of the resulting equations 
we give a concrete example by solving these equations numerically for the specific case of 
M-theory on AdS 4 xS 7 . 

Another motivation to consider a backreacted thermal gas on AdS^+ixS^ background, 
thus generalizing the work of (32], is in order to clarify the effect of the KK modes on the 
S fc sphere. In the previous sections of the paper we were working in the d + 1 dimensional 
reduced theory on AdS^+i, so a natural question was how much do the other KK modes on 
the sphere modify our results. To answer this question qualitatively it is perhaps easier to 
go directly to the full d + k + 1-dimensional theory where there is no issue about the KK 
reduction. As we will see the qualitative features of the solutions are similar to the ones we 
have been discussing. 

8.1 The equations of motion 

We start by considering a general AdS^+i x S fc compactification. We take the action in d+k + 1 
dimensions to be 



of Chern-Simons terms will not be relevant for our static spherically symmetric solutions. The 
equations of motion are 




where for the (d+l)-form we have defined F 2 



(d+iy. /w— w+i 



pv*i-Md+i_ The possible presence 





(8.1) 
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To fix the asymptotic form of the metric we first look for an AdS^ + ixS fc solution, without 
the radiation fluid, of the form 



ds 2 = fdsi dSd+i + R 2 dQ 2 k 

where ds\ dSd+i an< ^ d£l 2 k are the metrics of the unit radius AdS^+i and S k respectively. We 
take the form F to be proportional to the volume form of AdS^+i- We introduce the effective 
d + 1 dimensional cosmological constant A as 

F 2 (d + k-1) 



4 (d— 1)(A: — 1) 
and then from the equations of motion we find that 



A (8.2) 



d(d-l) p2 _ {d-l){k-lf 



2A 2dA 

In the presence of the radiation fluid the d + k + 1 dimensional equations of motion ( 18. II) 
are modified as 

Rfj,v ~ ~^9llvR — 2-jjF^ij,2...fid + i^v 2 fld+1 ~ 9iiv + ^GT^ 
where the stress energy tensor of the fluid in d + k + 1 dimensions is 

with the equation of state for radiation 

v=dTk? 

The energy density is related to the temperature in d + k + 1 dimensions by p = crT d+k+l 
where the constant a is proportional to total number of bosonic degrees of freedom of the 
theory. The velocity satisfies u fl u fl = — 1 with respect to the metric 'g iiv . 

We now focus on static configurations which also respect the SOik + 1) symmetry of the 
S h . The size of S k will not be constant in the radial direction, therefore we consider the 
following ansatz for the metric in d + k + 1 dimensions 

ds 2 = e-^'VM + R 2 e 2 ^dQ 2 k 

where g^ u is the d + 1 dimensional metric and the form F is still proportional to the volume 
form of the d + 1 dimensional space. To proceed we will write the equations in terms of the 
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d + 1 dimensional metric g and the scalar cf). The d + 1 dimensional Newton's constant G is 
defined by the relation 

1 R k V k 



16ttG 167TG 

where is the volume of the unit S k . The velocity is normalized to square to —1 with 
respect to the metric g^, so we introduce 



which satisfies g fiu u^u u = — 1. We also define effective d + 1 dimensional pressure and density 

p = (R'V^e-^p, p = {B!'V k )- 1 e-&+p 

After some algebra we find that ( 18. II) can be written as d + 1 dimensional equations in the 
form 

1 H<i + k — 1) ( 1 \ 

□0 = V'(4>) - 8nGp 

with 

a e d ~ lV e d - x v \ d — 1 



*W = -At— r -i-r- - -7^-— r - A- 



&d d + k-lj k(d + k-l) 
and 

8.2 The graviton star solution 

Now we consider a spherically symmetric ansatz for the d+1 dimensional metric 

ds * = _ e 2x{r)-2 P[ r) dt 2 + ^(r) ^2 + 

with -u t = e x( - r - ) ~ /3 ^ r - ) and <f)(r) , p(r) , p(r) functions of r only. The scalar equation becomes 
d 2 r <P + ( X ' - 2/3' + — ^ d r 4> = e 2/ V(0) ~ e^MGp 



The tt equation gives 

d ~ 1 ((d - 2)(e 2 ^ - 1) + 2r/3') e" 2 ^^-^ = ^ k ~ 1] (\{d r <j>fe^ + V{4>) ) +8nGp 



2r 2 u /v 7 ^ ' 2£ 2 (d-1) \2 
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Defining 

-2/3 



r 2 _ M(r) 
e -^ = i + __c 



£2 r rf_2 



with C d+1 = pjjj^ we find 



M ' = W "S.G (rf-1) V2 
Adding the and rr equations we find 

, _ r ( k(d + k- 1) 1 2 2/J 



Finally from the conservation equation of the total T Mi , one can show that the effective d+1 
dimensional density and pressure are 

p(r) = J^_ e -^^<l>n+{d+k+i)x a e M ^P 11 ^(r)-(d+fc+i)(x(0-/3W) ) p ( r ) = — *— p ( r ) ( 8 .4) 
R^Vk d ~\~ k 

where po is the d + k + 1 dimensional density at the center. 

These equations can be used to find solutions numerically. We have to specify the d + k + 1 
dimensional density po a t the center and the value <po of the scalar field and then integrate 
outwards. The value Xo = x(0) & t the center can be determined as follows: we notice that the 
differential equations are invariant under a constant shift of x which corresponds to an overall 
rescaling of the time coordinate. By demanding that at infinity the metric is written in the 
standard parametrization of AdS space we have to impose x{°°) = 0- Using the invariance 
under constant shifts of x this can always be achieved by an appropriate choice of xo- So 
Xo is not really an independent parameter that we can tune. At this stage it seems that 
the independent initial data are given by the two numbers (po,4>o)- However we also have 
another boundary condition that we have to satisfy. At large r the scalar field has to go to the 
minimum of the potential = 0, otherwise we would not have an asymptotically AdS space. 
This shows that the pairs of initial data (po, (po) are actually not independent. For every po 
there is a value 0o such that after evolving these initial data we get an asymptotically AdS 
space. Solving this "shooting problem" numerically we can determine a one-parameter family 
of graviton-star solutions parametrized by the density at the center p . 

In figure 1261 we show the mass of the star as a function of po- We notice that the qualitative 
features are similar to those of a radiation star [32J in AdS4. There is a critical central density 
p c beyond which the mass of the star starts to decrease signaling a Chandrasekhar-type 
instability. Solutions with p > p c are unstable under radial perturbations and presumably 
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Figure 26: Mass of a thermal star in AcISa X S dbS cl function of the 11-dimensional energy 
density po at the center. 



collapse. The natural endpoint of this collapse is an eleven dimensional Schwarzschild black 
hole in thermal equilibrium with its Hawking radiation^. It would be interesting to study 
these thermal stars in more detail and their possible relevance for finite temperature gauge 
theories. 



9 DISCUSSIONS 



In this paper we studied the holographic interpretation of degenerate fermionic stars in anti 
de Sitter space. Clearly our work leaves many important questions unanswered and is only 
a first step in a potentially interesting direction. Let us now mention some possible further 
directions. 

It may be interesting to study the dynamical process of collapse towards a black hole 
by computing numerically the time-dependent solution in the case that the mass of the star 
exceeds the Chandrasekhar bound. In principle this bulk computation is straightforward, 
though perhaps tedious, and it might be useful for answering certain questions about the dual 
deconfinement transition in the boundary CFT. Along these lines it might also be interesting 
to explore possible connections with critical scaling during the collapse. The Choptuik scaling 
has been discussed in the context of AdS/CFT in recent works [HHH]. The analogue of 



Choptuik scaling has been observed in the (driven) collapse of neutron stars [45, 



15 



SO 



it would presumably be relevant for the collapse of our star and might have an interesting 
interpretation in the dual gauge theory. 

Another interesting direction would be the computation of boundary correlation functions 
in the presence of the star. It would be useful to understand how the bulk Fermi surface 



44 We are always working in the microcanonical ensemble so the total energy is conserved. 

45 We would like to thank L. Alvarez-Gaume for bringing this to our attention and for useful comments. 
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is visible in terms of the boundary correlators and compare with works [T7T420] where the 
fermionic correlators are evaluated in the background of a black hole. Let us notice that our 
star (at least the one without charge) can only exist in global AdS coordinates and not in the 
Poincare patch. In terms of the boundary theory the star would describe a low-temperature 
sector of the theory, at temperatures below the Hawking-Page transition and the black hole 
formation. So in a sense it would be the analogue of the low-temperature phase of the systems 
studied in the aforementioned works, at low temperature/chemical potential. 

It would be useful to further clarify issues related to the validity of our approximations. 
While we tried to address some of these issues, a more carefully analysis would be desirable. 
In particular it it would be interesting to settle conclusively whether a degenerate fermionic 
star can be reliably realized in a known AdS / CFT duality. 

Finally, it is clear that the boundary CFT treatment in our paper was on a rather basic 
level. We hope that more can be done in the direction of analyzing the interacting gas of 
fermions directly on the CFT side, perhaps under suitable simplifying assumptions about the 
spectrum of the theory and the strength of interactions, extending the basic estimates that we 
sketched in this paper. Furthermore it would be interesting to understand how to translate 
the dynamical evolution of the collapsing star (after the relevant solution has been calculated 
in the bulk) into useful information about the thermalization and deconfinement of the gas 
of glueballs in the gauge theory side. We leave these questions for future work. 
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